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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


Amn 

solar  insolation  absorbed  in  the  layer  (m,n) 

a  (m,n) 

Manabe-Moller  absorptivity  function 

ABA 

absorptivity  of  the  troposphere 

ABG 

fractional  absorptivity  of  solar  insolation  by 
earth* s  surface 

ALB 

earth-atmosphere  system  albedo 

ATRAN 

transmissivity  of  the  troposphere 

Stefan-Boltzmann  blackbody  flux  at  T 

.K 

BALB 

24-hour  averaged  radiational  balance  at  earth's 
surface 

BALk1k2 

24-hour  averaged  radiational  balance  for  layer 
(kl'k2} 

BALT 

24-hour  averaged  radiational  balance  at  tropopause 

C 

carbon  dioxide  layer  absorber  mass 

cal  cm  min  calories  per  centimeter  squared  per  minute 


CL 

total  opaque  cloud  cover 

CLi 

fractional  cloud  amount  for  layer: 

I  =  1  in  600  to  400  mb; 

I  =  2  in  900  to  800  mb 

E 

East  longitude 

e 

X 

vapor  pressure  at  top  of  constant  flux  layer 

F  (A) 

solar  insolation  subject  to  water  vapor  absorption 
only 

F(2) 

effective  solar  insolation  at  tropopause 

FADJ 

total  incoming  insolation  at  top  of  atmosphere 

Fklk2 

net  infrared  flux  divergence  in  layer  (k^fk2) 
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V 

net  infrared  flux  at  level  k 

FNWC 

Fleet  Numerical  Weather  Central 

F  (S) 

solar  insolation  subject  to  Rayleigh  scattering  only 

ft 

net  IR  flux  to  space 

g 

2 

gravity  =  9.8067  m  sec 

f 

multiplicative  factor  for  tuning  cloud  reflectances 

h 

hour  angle 

H 

height  of  homogeneous  atmosphere;  24-hour  averaged 
hour  angle 

HL 

population  of  gridpoints  north  of  25N 

I 

abscissa  grid  location 

IA10(m,n) 

solar  insolation  absorbed  at  surface  with  cloud 
condition  (m,n) 

IS10(m,n) 

solar  insolation  at  surface  subject  to  Rayleigh 
scatter  with  cloud  condition  (m,n) 

J 

ordinate  grid  location 

k 

pressure  level  used  in  this  study  equal  to  10a 

i  •  -1 

ly  min 

langleys  per  minute 

N 

North  latitude 

P 

pressure 

pk 

pressure  in  millibars  (mb)  at  level  k 

qk 

mixing  ratio  at  level  k 

Qklk2 

24-hour  averaged  solar  warming  in  layer 

QAVE 

24-hour  averaged  insolation  at  the  tropopause 

solar  net  insolation  at  level  k  =  0 

r 

Bowen  ratio;  actual  earth-sun  distance 

r 

m 

mean  earth-sun  distance 
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R 

net  radiation  balance  at  the  surface 

R 

a 

mean  radiative  cooling  rate  in  troposphere 

Rd 

universal  gas  constant 

REF 

total  insolation  reflected  back  to  space 

REFA 

F(A)  insolation  reflected  back  to  space 

REFS 

F (S)  insolation  reflected  back  to  space 

RH 

relative  humidity 

Rt 

mean  radiative  energy  gain  (loss)  rate  at  ocean- 
troposphere  system 

rn 

radiative  net  flux  at  level  k  =  0 

s 

South  latitude;  effective  solar  constant 

s 

o 

heat  storage  in  oceanic  water  mass 

S 

a 

heat  storage  term  for  the  troposphere 

Study  A 

Spaeth's  Thesis;  using  winter  data  (see  references) 

Study  B 

Meyers'  Thesis;  using  spring  data  (see  references) 

Study  C 

Beahan's  Thesis;  using  summer  data  (see  references) 

Study  D 

Warner's  Thesis;  using  autumn  data  (see  references) 

Tk 

temperature  at  level  k 

TR 

population  of  gridpoints  20S-to-25N  inclusive 

TRAN 

total  insolation  incident  at  the  earth's  surface 

T 

X 

temperature  at  the  top  of  constant  flux  layer 

U 

water-vapor  layer  absorber  mass 

W 

West  latitude 

W  (m,n) 

cloud  fractional  weight  for  cloud  condition  (m,n) 

Z 

Zenith  angle 

a  (G) 

surface  albedo 
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a(R) 

6 

A (ALB) 
e 

WC 


A 

7T 

P 

a 


Rayleigh  clear  sky  albedo 

solar  declination  angle 

difference  between  ALBMOD  and  ALBRAS 

emissivity  due  to  water  and  carbon  dioxide  absorber 
mass  at  indicated  layer 

potential  temperature  at  level  k 

longitude 

surface  pressure?  pi  =  3.1416 
density 

sigma  pressure  level  used  by  FNWC,  normalized 
to  surface  pressure 

latitude 
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I.  INTRODUCTION 


This  thesis  is  a  refinement  of  previous  radiational  models  described 
by  (A)  Spaeth  (1975),  (B)  Meyers  (1975),  (C)  Beahan  (1975)  and  (D) 

Warner  (1974)  for  use  in  the  Fleet  Numerical  Weather  Central  (FNWC)  pre¬ 
diction  system.  This  study  has  as  a  primary  objective,  the  comprehensive 
re-examination  of  the  radiational  physics  in  layers  comprising  the  ocean- 
atmosphere  system.  FNWC  atmospheric  soundings  defined  at  constant  pres¬ 
sure  levels  for  the  four  mid-seasonal  dates  of  the  "data  year"  1973-74 
as  previously  examined  in  studies  (A,B,C,D)  were  utilized  in  this  study. 

The  application  of  the  FNWC  radiative  model  may  be  made  to  any  scale 
of  analysis  for  which  there  is  adequate  resolution  of  the  temperature 
and  moisture  data  in  the  vertical.  In  the  horizontal,  the  reliability 
of  the  data  used  here  is  consistent  with  that  of  the  FNWC  interpolation 
to  gridpoints  in  the  analysis  procedure.  Temperature  and  dew-point  data 
in  radiative  soundings  are  typically  reported  to  the  nearest  tenth  of  a 
degree.  The  radiative  computations  made  here  are  applied  to  FNWC  grid- 
points  and  are  designed  to  make  a  one-hour  forward-time  step  at  gridpoints 
in  the  FNWC  primitive  equation  forecast  model,  with  special  adaptions  to 
0-levels . 

The  specification  of  cloud  amounts  in  two  designated  layers,  one  at 
a  mid-level  and  the  other  at  a  low-level,  has  an  important  influence  on 
the  radiative-model  dispositions,  both  in  the  short-  and  long -wave 
spectral  regions.  Initially  in  (A,B,C,D),  the  specification  of  the 
fractional  amounts  of  CL^  and  CL^  had  been  based  on  large-scale  formula¬ 
tions  developed  by  Smagorinsky  (1960) ,  but  during  the  course  of  these 
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studies  it  was  found  more  realistic  to  modify  initial  CL-values  to 
CL*  =  2/3  CL.  This  reduction  in  cloud  amount  was  also  used  here  as  it 
had  been  chosen  to  prevent  CL*  from  exceeding  unity  and  to  afford 
better  agreement  with  satellite  climatology,  such  as  planetary  albedo  , 
for  compatible  data  periods . 

With  the  reduced  cloud  coverages  CL^ '  and  CL^ 1  from  the  earlier 
studies  (A,B,C,D),  it  was  possible  to  obtain  reasonably  close  agreement 
in  the  computed  terrestrial  net  flux  at  the  top  of  the  atmosphere  and 
that  observed  for  comparable  NIMBUS  III  subsatellite  points  and  data 
periods  (Raschke,  Von  der  Haar,  Bandeen  and  Pasternak,  1973).  However 
even  with  the  reduced  cloud  coverages  CL^ 1  and  the  computed 

planetary  albedo  remained  generally  excessive,  particularly  in  tropical 
latitudes.  Hence  it  became  a  major  objective  of  this  particular  study 
to  modify  empirically  the  reflective  capability  of  the  cloud  layers. 

This  was  done  by  defining  a  general  factor  f  so  that  the  initial  choices 
(after  Rodgers,  1967)  of  cloud  reflectances  R  were  modified  to  R*  where 

R*  =  fR  o 

Systematic  substitution  of  the  cloud  reflectances  R,  wherever  they 
entered  the  solar  disposition  equations,  by  R*  then  led  to  a  relation¬ 
ship  between  the  global  albedo  and  f.  Utilization  of  the  least  squares 
technique  to  minimize  the  differences  between  satellite  and  model  albedos 
over  a  geographic  sample  of  points  led  to  best-fit  value  of  the  "tuning- 
factor"  f.  Separate  values  of  f  were  deduced  by  least  squares  for  each 
season,  and  subselections  were  deduced  for  the  tropical  and  extra tropical 
areas ,  respectively . 
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The  modified  solar  cloud-reflectances  improved  the  agreement 
between  satellite  and  radiative  model  albedos  in  both  geographic 
areas  insofar  as  net  incoming  insolation  was  concerned.  The  terres¬ 
trial  net  flux  at  the  top  of  the  atmosphere  was  unaffected  by  the 
choice  of  f ,  while  the  use  of  CL^ '  and  CL^1  as  specified  gave  good 
agreement  with  satellite  terrestrial  net  flux  data  over  the  geo¬ 
graphic  range  and  the  time-scales  concerned. 
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II.  DATA  PREPARATION 


A.  DATA  FIELDS 

1.  General  Considerations 

The  initial  temperature  and  humidity  data  used  in  this  study 
were  arranged  in  the  form  of  soundings  taken  along  four  oceanic  meri¬ 
dians  (Fig.  1)  of  the  Fleet  Numerical  Weather  Central  (FNWC)  Northern 
Hemisphere  mid-seasonal  analyses  for  16  October  1973,  16  January  1974, 
16  April  1974  and  16  July  1974.  Oceanic  locations  for  these  computa¬ 
tions  were  chosen  because: 

p 

(a)  constant  a-surfaces  (where  0  =  — )  of  the  FNWC  primitive 
equation  system  are  nearly  identical  constant  pressure  levels. 

(b)  The  maritime-area  soundings  are  more  likely  to  be  system¬ 
atically  representative  of  the  set  of  zonally-distributed  gridpoints 
than  over  land. 

The  three  meridians  selected  over  the  Pacific  Ocean  were  lo¬ 
cated  at  125W  (25  soundings),  170W  (25  soundings)  and  145E  (17  sound¬ 
ings)  .  The  Atlantic  Ocean  meridian  was  located  at  35W  (26  soundings) . 
This  method  of  selecting  "soundings"  along  the  indicated  meridians 
of  the  FNWC  polar  stereographic  map  made  it  unnecessary  to  employ  spa¬ 
tial  interpolation  between  original  data  gridpoints  along  the  meri¬ 
dians.  Data  along  line  3  in  the  Pacific  was  not  extended  southward  of 
gridpoint  (9,55)  because  they  fell  over  land  masses  (New  Guinea  and 
Northern  Australia)  where  the  surface  temperatures  and  other  sounding 
features  were  unrepresentative  of  the  oceanic  values. 
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'Figure  1.  FNWC  polar  stereographic  grid  and  meridians  (lines  1,  2, 

3 /  and  4)  selected  for  study.  The  longitudes  A  are  shown 
for  each  meridian  as  well  as  the  extent  considered  of  each 
'  meridian. 
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2.  Data  Treatment 


a.  Original  Soundings  and  Modifications 

The  gridpoint  soundings  were  taken  from  the  original  FNWC 
63-by-63  Northern  Hemisphere  analyses  of  T(p)  and  of  T  -  T^  (the  dew 
point  depression)  at  standard  pressure  levels  up  to  and  including 
p  ^  100  mb.  Examples  of  such  original  soundings  were  shown  as  Table 
1(a)  in  Meyers  (1975,  p.  24).  Subsequently  each  original  FNWC  sound¬ 
ing  was  transformed  — in  previous  studies  (A,  B,  C,  D)  of  this  series 
— into  what  has  been  termed  the  radiative  sounding  having  the  format 
shown  in  Table  I.  The  data  levels  of  the  radiative  sounding  contains 
essentially  the  five  FNWC  predictive  O-levels  (dotted  levels  in  Fig.  2) . 

At  each  gridpoint  selected,  the  original  FNWC  humidity 
soundings  were  given  in  the  form  of  five  dew  point  depressions  over 
the  analyses  levels  from  925  mb  to  400  mb.  At  the  surface  (level 
k  =  10)  the  standard  instrument  level  vapor  pressure,  ea^r>  was  trans¬ 
formed  into  the  surface  mixing  ratio,  q^,  by  means  of 

q_n  =  621.97  (e  .  /1000)  .  (2-1) 

lu  ai  r 

To  obtain  radiative  soundings  as  in  Table  I,  it  is  neces¬ 
sary  to  have  water  vapor  and  CC>2  absorber  masses  at  certain  required 
k-level  boundaries  (Fig0  2) .  All  radiative  soundings  in  this  study 
start  at  sea  level  with  the  approximation  of  surface  pressure 
7T  =  1000  mb.  Therefore,  the  eleven  k-levels  correspond  closely  to  the 
FNWC  levels  P  =  1000.,  900.,  800.,  ....,  200.,  100.,  0.0  mb  and  in 
turn  to  =  1.0,  0.9,  ....,  0.1,  0.0. 
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TABLE  I.  Example  of  a  radiative  sounding  at  gridpoint  (1,1)  for 
16  April  1974  with  mixing  ratio  listed  at  odd  k-levels  (Fig.  2) . 
Additionally,  water-vapor  and  CO-  absorber  masses  are  also  listed 
as  these  parameters  have  been  modeled  in  the  radiative  theory  pre¬ 
sented  in  this  study. 


Pressure 

(mb) 

Temp 

(°C) 

Mixing 

Ratio 

(gAg) 

Absorber 
Water  Vapor 

(gm/cra^) 

Masses 

CO  2 

(cm/ cm 

1000 

25.60 

17.10 

900 

16.09 

12.08 

800 

11.21 

2.26 

45.53 

700 

5.00 

6.21 

600 

-2.10 

3.23 

83.53 

500 

-11.50 

2.57 

400 

-23.20 

1.35 

3.55 

113.35 

300 

-38.00 

0.65 

200 

-56.60 

3.61 

133.99 

100 

-80.70 

0.03 

3.61 

138.67 

0 

-80.70 

3.61 

143.35 
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a=1.0 


Figure  2.  Five-layer  radiative  sounding  used  in  this  study.  Levels 
are  identified  by  their  values  on  the  k-scale,  while  layers  are  identi¬ 
fied  by  their  level  boundary  indices  in  parentheses,  e.g.  (8,10), 
Pressure-scaled  water  vapor  and  CO^  mass  increments  Au  and  AC,  respec¬ 
tively  are  integrated  relative  to  the  surface  and  the  resulting  U  and  C 
are  carried  at  even  levels.  The  temperature  T  is  retained  at  all  levels. 
Amounts  of  clouds  CLj  and  CL2  in  the  layers  shown  have  been  para¬ 
meterized  for  consideration  of  their  radiative  effects. 
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b.  Radiative  Temperature  Profiles 

The  gridpoint  temperatures  were  listed  at  each  mandatory 
level  of  Table  I  between  1000.,...,  100  mb  (i.e.,  between  k  =  10,..., 
1) .  The  temperature  was  assumed  to  be  isothermal  from  100  mb  to 
0.0  mb.  The  temperature  at  level  k  =  10  was  set  equal  to  the  FNWC 
listed  sea-surface  temperature.  The  radiative  sounding  temperatures 
for  the  remaining  k-levels  were  obtained  from  either  their  correspond¬ 
ing  listed  temperature-level  or  by  a  three-point  Lagrangian  interpola¬ 
tion  scheme  [Eq.  (2-1),  Spaeth,  1975],  to  level  k  when  the  listed 
FNWC  temperature  profile  did  not  include  the  value  T^. 

c.  Radiative  Moisture  Profiles 

Similarly,  the  moisture  profiles  of  Table  I  have  been  ob¬ 
tained  by  an  interpolative  procedure  over  the  original-level  FNWC  mix¬ 
ing  ratios  to  those  required  at  k-levels  in  a  manner  analogous  to  that 
discussed  by  Spaeth  (1975,  pp.  29-31) . 

d.  Pressure-Scaled  Absorber  Masses 

The  pressure-scaled  water  vapor  and  the  carbon-dioxide 
scaled  absorber  masses  were  calculated  for  the  six  even  numbered  k- 
levels  (Fig.  2)  using  Eqs.  (2-8,  2-9,  2-10  and  2-11)  respectively  as 
outlined  by  Spaeth  (1975,  pp.  31-33) .  These  equations  use  essentially 
the  mixing  ratios  of  water  vapor  and  of  CC>2  at  odd  k-levels. 

B.  CLOUD  PARAMETERIZATION 

The  relative  humidities  (RH)  at  levels  k  =  5  and  k  =  9  are  used 
in  the  calculations  of  the  fractional  cloud  covers  CL^  and  CL^  in 
layers  (4,6)  and  (8,9)  respectively  (Fig.  2).  The  equations  for 
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parameterization  of  the  two  fractional  cloud  amounts  are  as  follows: 


CL  =  2/3  [2.0  (RH ( 5) )  -  0.7] 


( 2-2a) 


CL2  =  2/3  [3.33 (RH (9) )  -  2.0] 


(2-2b) 


The  bracketed  part  of  the  equations  (after  Smagorinsky,  1960)  were 
reduced  by  the  1 /3  factor  in  an  attempt  to  tune  cloud  amounts  to  obtain 
albedo  values  in  closer  agreement  with  the  recent  satellite  radiational 
climatology  of  Raschke  et  al.,  (1973).  The  "2/3  CL"  parameterization 
set  forth  in  Eq.  (2-2)  is  used  here  for  estimating  large  scale  radia¬ 
tional  effects  only.  Thus  small-scale  convective  activity  and  ci  priori 
climatological  effects  were  not  considered  in  specifying  the  form  of 
Eq.  (2-2)  .  Tuning  of  the  model-albedo  values  of  this  study  by  varying 
cloud  reflectance  coefficients  with  respect  to  season  and  latitude  will 
be  discussed  in  Section  IV. 

C.  CLOUD-AREA  COVERAGES 

Fractional  cloud  amounts,  CL^  and  CL^,  were  computed  at  each  grid- 
point  for  levels  k  =  5  and  k  =  9  by  Eqs .  (2-2a,b) .  In  addition  the 
gridpoint  area  may  be  thought  of  as  broken  into  random  fractional  seg¬ 
ments  of  size 


W (0,0)  =  (1-CL1)  (1-CL2) 


(2-3a) 


wherein  there  is  a  combination  of  clear-over-clear  segments  in  the 
layers.  Similarly,  the  gridpoint  area  has  the  fractional  area  of  cloud 
coverage 


W(l,l)  =  CL1  *CL2 


(2- 3b) 
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of  an  upper  overcast  amount  overlying  a  lower  overcast.  Likewise  the 
area-combinations  of  overcast  over  clear,  and  clear  over  overcast 
areas  in  two  layers,  may  be  visualized  as  occurring  with  the  weights 

W (1 , 0)  =  CL  *(1-CL2)  (2-3c) 

and 

W(0,1)  =  (1-CL  )  *CL2  (2-3d) 

respectively . 

For  radiational  computations  it  was  useful  to  carry  the  relative 
weights  or  fractions  of  the  gridpoint  area  exposed  to  the  specified 
cloud- layer  combinations.  Henceforth,  the  symbols  denoted  by  W(0,0), 

W  (1 , 1) ,  W(1,0)  or  W(0,1)  indicate  the  fractionally  overcast  (1)  or 
clear  (0)  cloud-area  combinations  in  the  indicated  layers  (Fig.  2) , 
with  the  first  index  1  or  0  referring  to  layer  CLn ,  k  =  5,  and  the 
second  to  CL^,  k  =  9. 

The  usefulness  of  the  cloud-area  weighting  device  will  be  clarified 
in  Sections  III  and  IV,  where  the  procedures  for  the  terrestrial  and 
solar  radiational  computations  are  discussed  and  the  results  are 
summarized  over  the  set  of  soundings. 

A  measure  of  the  effective  cloud-cover  area  which  has  been  found 
useful  in  previous  radiational  studies  has  been  the  total  opaque  cloud 
cover,  CL,  referring  to  the  amount  of  thick  cloud  cover  overhead  re¬ 
gardless  of  the  level.  For  the  cloud  model  presented  here  CL  may  be 
expressed  as 

CL  =  CL1  +  CL2  -  CL  *CL2  .  (2-4) 
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III.  TERRESTRIAL  RADIATION 


A.  THEORETICAL  AND  EMPIRICAL  BASIS 

Empirical  formulas  were  developed  by  Sasamori  (1968)  for  flux 
emissivities  in  the  atmosphere  associated  with  computations  for  the 
radiative  balance  requirements  of  the  NCAR  General  Circulation  Model. 
Sasamori  derived  the  empirical  emissivity  formulas  for  v/ater  vapor 
and  CO^  by  comparison  with  the  theoretical  values  built  into  the 
Yamamoto  Radiation  Chart  (1952) .  The  Yamamoto  chart  has  proved  to 
be  quite  accurate  for  numerical  checks  of  the  Sasamori  emissivities. 
This  chart  was  also  used  in  the  previous  studies  (A,  B ,  C,  D)  as  a 
systematic  guide  for  integration  of  the  radiative  transfer  formulas 
developed  by  Martin  (1972,  1975) ,  who  adapted  the  Sasamori  emissiv¬ 
ity  formulas  to  the  particular  layers  of  interest  in  the  gridpoint 
computations  of  the  FNWC  primitive  equation  model  (Fig.  2) . 

The  essential  long-wave  (IR)  net-flux  parameters  required  for 

use  in  this  study  are  the  following: 

* 

F  =  IR  net  flux  at  earth,  k  =  10 

* 

F  =  IR  net  flux  at  level  k  =  8 

o 

* 

F  =  IR  net  flux  at  level  k  =  6 
6 

★ 

F.  =  IR  net  flux  at  level  k  =  4 
4 

* 

F^  =  IR  net  flux  at  level  k  =  2 
In  addition  the  IR  net-flux  divergence  coolings  to  be  computed  at 
each  gridpoint  are 

F810  =  IR  net-flux  divergence  in  the  layer  (8,10) 

F68  =  IR  net-flux  divergence  in  the  layer  (6,8) 
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F46  =  IR  net-flux  divergence  in  the  layer  (4,6) 

F24  =  IR  net-flux  divergence  in  the  layer  (2,4) . 

In  the  Radiation  Balance  Studies  in  the  series  A,  B,  C,  D  only 
F610  and  F26  were  computed  because  time  restraints  in  the  present 
FNWC  operational  heating  package  have  prevented  the  use  of  greater 
resolution  in  the  vertical.  Here,  four  flux  divergences  are  computed 
in  order  to  examine  more  closely  the  variability  of  the  flux  diver¬ 
gences  over  the  layer  thicknesses  reduced  to  approximately  200  mb 

each.  To  compute  these  four  flux  divergences  it  was  necessary  to 

*  * 

utilize  additional  formulas  for  F  and  F  as  developed  by  Martin 

o  4 

(1975) . 

In  order  to  make  IR  net-flux  calculations  along  the  path  of  in¬ 
tegration,  there  must  be  a  physically  sound  representation  of  the 

emissivity  (£  )  as  a  function  of  both  water  vapor  and  C0_  absorber 

wc  2 

masses  in  layers  along  the  sounding.  For  a  complete  discussion  of 
the  emissivity  formulas  used  in  the  quadrature  scheme,  refer  to 
Spaeth* s  Appendix  A  (1975) . 

B.  NET  FLUX  FORMULATIONS 

The  radiative  sounding  as  depicted  in  Table  I  was  computed  as 

the  combination  of  parameters  U(k,10),  C(k,10)  and  T  for  each  re- 

K 

quired  level,  k  =  10,  8,  ...  1,0.  Cloud  parameters  CL^  and  CL2  were 

also  computed  by  Eq.  (2-2)  at  each  gridpoint  and  in  general  are  both 

non-zero.  The  grid  area  was  then  considered  to  be  composed  of  areal 

fractions  (weights)  defined  in  Eqs.  (2-3a,b,c,d)  and  denoted  by  the 

symbols  W(0,0) ,  W(l,l),  W(1,0),  W(0,1). 

* 

The  composite  net  flux  F^Q  (CL^,  CL^)  at  level  k  =  10  at  each 
gridpoint  is  then  constructed  by  using  the  appropriate  weight  factors 
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•k 

to  multiply  the  reference  net  flux  F^  computations  defined  for  the 
four  special  cloud-cover  cases  already  defined  in  Section  II. C: 

F10*(0'0)'  F10*(1'0)'  F10*(0'1)'  F10*(1'1) 

It  therefore  follows  that 


F10*(CL1'  CL2)  =  W(°'0)F10*(0'0)  +  W(1'°)F10*(1'0) 

k  k 

+  W(0,1)F10  (0,1)  +  W(1,1)F  (1,1).  (3-1) 

★ 

The  reference  net  fluxes  of  Eq.  (3-1)  are  associated  with  (1) 

clear  skies  in  both  layers,  (2)  overcast  in  the  upper  layer  only, 

(3)  overcast  in  the  lower  layer  only  and  (4)  overcast  in  both  layers 
respectively. 

Spaeth  (1975)  has  listed  these  reference  net  flux  formulations 
in  his  Eqs .  (3-6),  (3-7)  and  (3-8).  Using  the  definitions  of  W(0,0), 

W(1,0),  W(0,1)  and  W(l,l),  F^*  (CL^,  CL^)  Can  s^own  to  assume  the 
form 


F10*<CL1'  CL2>  =  [1-CL2  ^  ^Fwc  ^8' 1C1^  ^B10_B8^ 

+  (ewc(8,10)  +  Ewc(6,10))(B8-B6)]} 

+  (l-CLo  ) (1-CLh  ){b  -.5[(£  (6,10) 

+  £wc(4'10))(VB4)  +  (£wc(4'10) 

+  £wc(2,10))(B4-B2)  +  (Ewc(2,10) 

+  ewc(l,0))(B2-Bl)  +  Ewc((0,10)/T1)*B1]} 

+  CL9  {  (B  — B  )  [1.5c  (9,10)]}  .  (3-2) 

10  9  WC 

Here 


B 


k 


1.170403  x  10  7  T  4 
k 


(3-3) 
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is  the  Stefan-Boltzmann  blackbody  flux  in  langlies  per  day. 

Further,  £  (U,  ,C .  ,10)  is  the  combined  water-vapor  and  CCt  emissiv- 

wc  k  k  2 

ity  along  the  path  from  level  10  to  level  k.  This  emissivity  is  con¬ 
sidered  by  Sasamori  to  be  temperature  independent  for  T  >_  210K,  whereas 

£  represents  the  temperature  dependent  emissivity  applicable  for 
wc 

T  <  210K.  [See  pp.  136-137,  Spaeth  (1975);  Sasamori  (1968)]. 

* 


The  formulas  for  F  (CL  ,  CL  )  ,  k  =  2,4,6,  and  8  have  been  developed 

.K  _L  2 

by  Martin  (1975)  in  a  manner  analogous  to  the  derivation  of  the  weighted 
* 

F  .  The  results  are  reproduced  as  the  following  equations: 


F„  = 


[1 

-CL1 

] 

{B8 

-.5[e  (6 

wc 

,8) 

(V 

-v 

+ 

<£wc 

(6 

,8) 

+  £ 

wc 

(4, 

8)) 

(V 

-v 

+ 

(£ 

(4 

,8) 

+  £ 

(2, 

8)) 

(B.- 

-B„) 

wc 

WC 

4 

2 

+ 

(£ 

wc 

(2 

,8) 

+  £ 

WC 

(1, 

8)) 

(B2- 

"V  + 

v/c 


r  l 


+  CL±  [l-.5ewc(6,8)]  (BQ-B6)+  CLX  (l-CL^  [  1-  .  5^  (8, 10)]  (B^-Bq) 


+  (1-CL.) (1-CL0) [1-.5S  (8,10)] (Bn  -B_) 

1  2  wc  10  8 


(3-4) 


F6*  -  tl-CLj)  (B8-.5[euc(6,8)(B8-B6)  +  (4, 6)  (B^) 

+  <EWC<4'6>  +  £wc<2'6))  (VB2>  +  <E„c<2'6> 

+  £wc(l,6)) (B2-B1)  +  £wc((0,6) ,T1)*B1]} 

+  (1-CL. ) (1-CL.) {(Bin-BQ) [l.-.5(e  (6,8) 

+  £  (6,10))]}  +  CL.  { (B  -B  ) * 

WC  18  6 

[1-.5E  (6,8) ] }  +  CL  (1-CL-)  { (B  -B  ) * 

wc  1  2  10  8 

[1-.5(E  (6,8)  +  £  (6,10))]}  .  (3-5) 

wc  wc 
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f4*  =  f1-CLi]  Vt5[£wc(4'6)(VB4) 

+  (£wc(4'6)  +  £wc(4'8))(VV  +  £wc(2'4)(VB2) 

+  (£wc(2,4)  +£wc(1,4))(B2-B1)  +  ewc((0,4),T1)*B1]} 

+  CL.  {B  . 5 [£  (2 , 4)  (B  -B_)  (3-6) 

1  4  WC  4  z 

+  (£  (2,4)  +  e  (1,4))(B  -B  )  +  e  (  (0,4)  ,T  )  *B  ] } 

WC  WC  Z  1  WC  1  1 

+  (1-CL1)  (1-CL2)  {(B10-B8)  [l-.5(ewc(4,8)  +  ^ (4 , 10) ) ] } 

F2*  =  [1"CL1]  {B8"*5[£wc(2'4)(B4_B2)  +  (£wc(2'4) 

+  £wc(2'6))(VV  +  (£wc(2'6)  +  £wc(2'8))* 

(B  -B  )  +  E  (1,2) (B  -B  )  +  £  ( (0,2) ,T  )*B  ] }  (3-7) 

8  6  v/c  z  1  wc  1  1 

+  (l-CL  ) (1-CL  )  {(B.  -Br) [l-.5(Et  (2,8) 

1  z  10  o  WC 

+  £wc(2'10))]}  +  CL1  {B4_'5[£wc(2'4) (B4"B2) 

+  £  (1,2) (B  -B  )  +  £  ((0,2),T  )*B  ]}  . 

WC  Z  1  wc  1  1 

As  was  described  by  Spaeth  (1975,  Section  III.B.5.)  concerning  the 

use  of  the  composite  case,  Eq.  (3-2)  can  be  reduced  to  give  expressions 
* 

for  for  the  various  reference  cloud-cover  cases  (0,0),  (1,1),  (1,0) 

*  * 

and  (0,1).  The  resulting  schematics  in  the  case  of  F  and  F  are 

8  4 

depicted  as  the  unhatched  area  in  Figs.  3(a,b,c,d)  and  4(a,b,c,d), 

*  *  * 

respectively  below.  Similar  graphs  for  F^  ,  F^  and  F^  can  be  found 
on  pp.  41-43  of  Spaeth  (1975) ,  and  are  not  reproduced  here. 

* 

A  typical  gridpoint  listing  of  the  IR  net-flux  computations,  F 
(k  =  10, 8,... 2),  has  been  reproduced  in  Table  II  for  gridpoint  (1,1) 
based  upon  the  radiative  sounding  of  16  April  1974  (see  Table  I) .  The 
printout  procedure  involves  computation  of  the  reference  net-flux 
values 

Fk*(0,0),  Fk*(l,0),  Fk*(0,l),  Fk*(l,l) 
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0  0 
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Figure  3C  Terrestrial  net  flux  F8*  with  (a)  clear  skies  (case  (0,0));  (b)  high  overcast  only 
case  (1,0));  (c)  low  overcast  only  (case  (0,1));  (d)  both  high  and  low  overcast 
(case  (1,1)  )0  Unhatched  area  in  each  case  depicts  F  *«, 


0.0 
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Figure  4.  Terrestrial  net  flux  F  *  with  (a)  clear  skies  (case  ),0);  (b)  high  overcast  only 
(case  (1,0));  (c)  low  overcast  only  (case  (0,1));  (d)  both  high  and  low  overcast 
(case  (1,1)).  Unhatched  area  in  each  case  depicts  F  *„ 


from  each  of  Eqs.  (3-2),  (3-4),  (3-5),  (3-6)  and  (3-7)  for  level 

k  =  10,  8,... 2  respectively.  Then  Eq.  (3-1)  with  the  appropriate 

weight-factors  of  Eqs.  (2-3a,b,c,d)  has  been  utilized  to  derive  the 
* 

composite  F  (CL  ,  CL  )  values  that  are  listed  on  the  bottom  line  of 

jC  -L  ^ 

Table  II. 


TABLE  II.  A  sample  listed  of  IR  net-flux  computations,  weighting 
factors  and  composite  values,  F  * (CL  ,  CL2) ,  as  computed  for  grid- 
point  (1,1)  for  16  April  1974.  KNet  flux  values  in  ly  min'”1. 


Cloud  Case 

Weight 

(CLj^  CL2) 

W(CL1, 

CL2) 

* 

Fio 

* 

F8 

* 

F6 

* 

F 

4 

* 

F2 

(0,0) 

.1002 

.1664 

.2288 

.2649 

.3160 

.3487 

(1,0) 

.1497 

.0943 

.1331 

.0915 

.1706 

.2338 

(0,1) 

.3007 

.0512 

.1563 

.2341 

.2871 

.3202 

(1,1) 

.4494 

.0512 

.0586 

.0607 

.1706 

.2338 

* 

Fv  -composite  values 

.0692 

.1171 

.1379 

.2202 

.2713 

C.  TROPOSPHERIC  COOLING  BY  LAYERS 

* 

?  Fio 

COOLING 

For  each 

mid-seasonal  ■ 

day  listed  and  at  each  gridpoint, 

an  IR  net 

flux  computation  in 

the  format  of 

Table  : 

El  is  easily 

converted  into 

four  sets  of 

layer  cooling 

effects 

• 

* 

F810  =  F 

o 

* 

-  F 

10 

F68 

* 

=  F6 

* 

‘  F8 

F46 

* 

=  F4 

* 

~  F6 

F24 

* 

=  F2 

* 

-  F4 

(3-8) 

These  layer 

cooling 

rates 

(ly  min 

^  have 

then  been  collected 

in  meri- 

dional  cross-section  format  for  each  longitude  under  study  (see  Section 
V)  and  by  season. 
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The  overall  tropospheric  cooling  rate  by  IR  net  flux  is  then  given 
*  * 

simply  by  -  F^Q  at  each  gridpoint  for  the  date  of  the  radiational 

sounding.  The  tropospheric  cooling  rates  computed  are  then  identical 
to 

★  * 

F2  -  F  =  F24  +  F46  +  F68  +  F810  .  (3-9) 

★  * 

The  values  of  -  F^  so  deduced  are  discussed  on  both  a  seasonal 

and  a  zonally-averaged  basis  in  Section  V. 

It  will  suffice  to  discuss  here  the  zonally-averaged  values  of 

* 

F  (CL^,  CL^)  as  computed  by  the  long-wave  radiational  model  previously 

* 

presented  in  this  sectionQ  These  results,  listed  simply  as  F  in 
Table  III,  will  be  discussed  as  a  function  of  seasons,  latitudes  and 

* 

CL  (total  opaque  cloud  cover  given  by  Eq.  (2-4)).  The  listings  of  F^ 

in  Table  III  are  essentially  as  extracted  from  computations  in  the 

* 

format  of  Table  II  followed  by  meridional-averaging  of  F  across 

constant  latitude  lines.  Finally  the  zonally-averaged  annual  values 
* 

of  both  F^q  and  CL  have  been  computed  by  arithmetic-averaging  over 
the  four  mid-seasonal  results  at  each  five-degree  increment  of  latitude 
from  20S  to  65N  (cf.,  Eq.  (5-1)). 

* 

The  model-annual  values  of  F^  presented  in  Table  III  are  presented 

with  those  derived  from  Budyko  (1956) ,  which  in  turn  are  listed  in  the 

final  column  of  Table  III.  Corresponding  values  of  total  opaque  cloud 

cover,  CL,  for  the  Budyko  climatology  were  not  available  so  that  only 

★ 

a  general  comparison  of  the  two  annual  F  zonally-averaged  distribu¬ 
tions  is  possible. 
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Annual 

16  January  16  April  16  July  16  October  Model  Values  Budyko 
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* 

Table  III  depicts  the  zonally-distributed  values  of  and  of  CL 

at  the  earth's  surface.  In  a  seasonal  comparison  of  the  model-computed 

*  * 

values  it  is  clearly  shown  that  (CL)  i^s  ci  decreasing  function 

of  cloud  cover.  There  is  a  clear-cut  tendency  in  each  season  for  a 

★ 

maximum  value  of  F^Q  to  be  located  in  the  subtropics  (latitudes  15N- 

* 

25N) .  Also  there  is  evidence  of  a  high  latitude  (55N-60N)  minimum  F^ 
associated  with  a  concentration  of  maximum  cloud  cover  CL.  An  outstand¬ 
ing  variation  is  the  transition  in  the  Southern  Hemisphere  latitudes 
(20S  -  10S) ,  which  has  small  cloud  cover  in  local  winter  and  compari- 
tively  large  cloud  cover  during  the  other  three  data  periods.  This 

cloud-cover  variation  corresponds  in  general  to  the  ITCZ  behavior  in 

* 

these  latitudes  across  the  indicated  seasons;  so  that,  F^  is  a  maximum 

in  mid-July  and  a  relative  minimum  in  the  period  January- April . 

The  final  entry  in  Table  III,  namely  "Wt.  Avg." ,  is  the  Northern 

Hemisphere  mean,  cosine-weighted  with  respect  to  latitude  (cf.,  Eq.  (6-1)). 

* 

The  cosine-weighted  F  values  show  a  minimum  in  Northern  Hemisphere 
summer  with  no  clear-cut  differences  in  CL.  This  summer  minimum  is 
attributable  to  higher  downward  flux  with  increased  summer  vapor 
pressures . 

D.  OUTGOING  IR  NET  FLUX  TO  SPACE 

1 .  Parameterization  Formula  for  Top-of-Atmosphere  IR  Net  Flux 

In  the  four  previous  radiation  studies  (A,  B,  C,  D) ,  an  approxi¬ 
mation  to  the  IR  net  flux  to  space,  designated  as  FF2 ,  was  computed 

from  the  radiative  sounding  at  each  gridpoint.  FF2  was  essentially 

★ 

an  extrapolation  of  F^  to  the  top  of  the  atmosphere  obtained  by  delet¬ 
ing  the  downward  IR  flux  due  to  stratospheric  water-vapor  and  CO^.  A 
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more  precise  expression  for  the  net  IR  flux  to  space  was  introduced  here 
after  Martin  (1975) o  The  development  is  analogous  to  the  quadrature 


summations 


B=0 


F_  =  /  e  dB 

T  J  wc 

Bn  „ 


(3-10) 


for  the  various  reference-cloud  combinations  (0,0),  (1,0),  (0,1)  and  (1,1) 

of  Fig.  5(a,b,c,d)  respectively .  The  final  quadrature-formula  is  then 
obtained  as  the  weighted  net-flux  result  as  was  also  done  in  Eqs.  (3-4)  , 
(3-5) ,  (3-6) ,  (3-7)  and  is  listed  below. 


FT  =  [1~CL1] {Bq-,5 [ (^wc (Of  1)  +  ewc(0'2))  (B2-Bl) 

+  (£  (0,2)  +  £  (0,4) ) (B  -B  ) 

wc  wc  42 

+  (£  (0,4)  +  e  (0,6))  (B-B.)  +  (e  (0,6)  +  £  (0,8))  (B  -B  )]} 

wc  wc  6  4  wc  wc  8  6 

+  CL  {B-.5[(£  (0,1)  +  £  (0,2) ) (B-B  )  (3-11) 

1  4  wc  wc  2  1 

+  (e  (0,2)  +  £  (0,4)  )  (B  -B  ) ]} 

wc  wc  4  2 

+  (l-CL  )  (1-CL  )  {  [1- .  5  (E  (0, 8)  +  £  (0,10))]  [B-B  ]} 

1  2  wc  wc  10  o 


where  CL  ,  CL^  are  given  by  Eq.  (2-2)  .  Note  that  F^  of  Eq.  (3-11)  is 
representable  by  the  unhatched  areas  of  Figs.  5(a,b,c,d)  for  the  various 
reference-cloud  cases  and  that  F^  has  no  downward  IR  flux  corresponding 
to  the  level  k  =  0. 

2 .  Comparisons  of  F^  with  Mid-Seasonal  Satellite  Climatology 

Comparison  was  made  of  computed-model  values  of  F^  with  satellite 
measurements  of  total  long-wave  flux  to  space  (after  Raschke  et  al.,  1973) 
for  the  Nimbus  III  mid-seasonal  periods  most  nearly  comparable  to  that 
of  the  FNWC  data.  F^  for  each  mid-seasonal  date  was  computed  over  the 
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Figure  5.  Outgoing  terrestrial  net  flux  to  space  with  (a)  clear  skies  [case  (0,0)];  (b)  high  ove 

cast  only  [case  (1,0)];  (c)  low  overcast  only  (0,1);  (d)  both  high  and  low  overcast  [case 

(1,1)].  Unhatched  area  in  each  case  depicts  F  . 


four  meridians  considered.  These  Ft  values  were  then  averaged  across 
the  four  meridians  to  get  a  mean  zonal  distribution  of  the  type  shown 
in  Table  IV. 

Table  IV  shows  the  zonally-averaged  model  values  compared  with 
those  extracted  from  Raschke.  Raschke's  results  were  obtained  by 
averaging  across  the  same  oceanic  meridians  as  those  used  in  this  study. 
Again,  in  the  bottom  line  of  each  column  in  Table  IV  is  listed  the 
cosine-weighted  mean  of  each  set  of  column  values  for  the  Northern 
Hemisphere  only. 

The  zonally-averaged  values  computed  by  the  F^-model  are  very 
close  to  those  reported  by  Raschke,  especially  between  0-65N  and  in 
the  Northern  Hemispheric  means .  The  limitations  of  the  comparisons 
made  here  are  obvious,  when  it  is  recalled  that  between  latitudes  20S- 
5S  and  between  60N-65N  there  are  fewer  than  four  meridional  lines 
available  for  computing  the  listed  zonal  values  in  Table  IV.  For  all 
other  zonally-averaged  values,  four  meridional  lines  were  used  in  the 
averaging. 

The  close  comparison  between  model-values  of  F^  from  (3-11)  and  those 
essentially  derived  from  satellite  climatology  tend  to  support  the 
cloud  parameterization,  Eq.  (2-2) ,  insofar  as  IR  net  flux  is  concerned. 
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IV.  SOLAR  RADIATION 


A.  COMPOSITION  OF  SOLAR  INSOLATION 

At  the  top  of  the  atmosphere  (k=0)  this  study  assumed  a  solar  con¬ 
stant  of  2.00  ly  min  1  (Joseph,  1971).  Furthermore,  this  constant  was 
assumed  subject  to  a  four  percent  attenuation  above  the  tropopause  due 
to  ozone  and  oxygen.  Thus  the  effective  solar  constant  at  level  k=2 
in  this  study  is  1.92  ly  min  1. 

To  compute  the  effective  solar  insolation  at  the  tropopause  the 
following  formula  was  used 

F  (2)  =  S  [— — ]  ” 2  Cos  Z  (4-1) 

m 

where  S  is  the  effective  solar  constant  at  level  k=2  and 

r/r  =  ratio  of  the  actual  earth-sun  distance  to  the  mean 
m 

earth-sun  distance,  a  function  of  the  Julian  date 
Cos  Z=  cosine  of  the  zenith  angle,  a  function  of  the 
Julian  date  determined  by 

Cos  Z  =  Sin  <j)  Sin  6  +  Cos  <j)  Cos  6  Cos  h  .  (4-2) 

The  symbols  on  the  right  side  of  Eq.  (4-2)  are  defined  as  follows: 

(J)  =  latitude 
6  =  solar  declination 
h  =  hour  angle 
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The  Smithsonian  Meteorological  Tables  (List,  1958)  list  the  ratio  of 


r/r  and  the  solar  declination,  6,  for  the  mid-seasonal  dates  appli- 
m 

cable  to  this  study  and  reproduced  in  Table  V. 


TABLE  V.  Values  of  the  ratio  of  the  earth-sun  radius  vector,  r/r  , 
and  of  the  solar  declination  angle,  6,  used  in  this  study  for  year 


1973  -  1974. 

Date 

r/r 

m 

6 

16  January 

0.98372 

21 . 07917°S 

16  April 

1.00333 

8 . 48333°N 

16  July 

1.01644 

21 . 50000°N 

16  October 

0.99717 

8. 22500°S 

The  value  of  sin  (p  was  calculated  using  one  of  two  different  formulas, 
depending  on  the  data-line  used  for  the  computations,  in  terms  of  the 
FNWC  map  coordinates  (I,J)  as  in  Eq.  (4-3a,b) .  Conversely  for  these 
lines  one  may  solve  for  I  in  terms  of  sin  (p  as  in  (4-3c,d)  : 


Lines  1,3,4 


Line  2 


Lines  1,3,4 


Sin  (p 


Sin  (p 


973.752  -  2 (32-1) 2 

973.752  +  2 (32-1) 2 
973.752  -  (32-1) 2 

973.752  +  (32-1) 2 


I  =  32  -  22.065  [ 


Cos  (p 
1  +  Sin  (p 


(4-3a) 


(4-3b) 


(4-3c) 


Line  2 


31.205 

,25 

,25 

,38 


Cos  (p 
1  +  Sin  <J) 

for  Lines  1,2 
for  Line  3 
for  Line  4 


(4- 3d) 


Here  I  is  the  abscissa  distance  on  the  FNWC  grid  (Fig.  1)  and  varies 
by  line  as  described  in  Section  II.  The  soundings  for  lines  1,  2,  and 
3  were  all  taken  at  OOOOGMT  with  the  solar  noon  existing  at  the  180th 
meridian;  therefore  the  hour  angles  for  these  three  lines  were  55°, 

10° ,  and  35° ,  respectively.  For  line  4,  the  soundings  were  taken  12 
hours  earlier  with  solar  noon  at  the  Greenwich  meridian,  giving  an 
hour  angle  for  line  4  of  35°. 

A  very  simple  partition  of  solar  insolation  was  utilized  in  this 
study  after  Joseph  (1971) .  It  consisted  of  dividing  the  insolation 
F(2)  into  two  parts  at  level  k=2.  One  part  was  considered  to  include 
all  wavelengths  X  >  .  9  ym  where  absorption  by  water  vapor  and  carbon 
dioxide  bands  are  the  most  prevalent  attenuation  processes  in  clear  air. 
This  part  of  the  solar  spectrum  was  termed  the  F (A)  energy  and  con¬ 
sidered  subject  to  water-vapor  absorption  but  not  to  Rayleigh  scatter¬ 
ing.  For  shorter  wavelengths,  X  <  .9  ym,  absorption  of  the  solar  in¬ 
solation  energy  by  water  vapor  was  considered  negligible.  This  part  of 
the  solar  insolation  was  denoted  F(S)  suggestive  of  the  fact  that  it 
was  subject  only  to  Rayleigh  scattering  attenuation  in  clear  air.  The 
two  partitions  are  formulated  after  Joseph  (1971)  as  follows: 


F  (A) 

=  .349  F (2) 

(4-4) 

F  (S) 

=  .651  F (2)  . 

(4-5) 

In  this  study,  the  introduction  of  two  cloud  decks  produced  cloud-reflec¬ 
tivity  effects  upon  both  the  F (A)  and  F(S)  solar  energy  insolations. 
However,  in  the  clear  areas  around  any  gridpoint  only  the  absorption- 
attenuation  applies  to  the  F (A)  insolation,  while  only  Rayleigh  scatter¬ 
ing-attenuation  applies  to  the  F(S)  insolation. 
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B.  DISPOSITION  OF  F(S)  INSOLATION 


In  the  disposition  of  the  F(S)  insolation,  Joseph  (1971)  determined 
that  Rayleigh  scattering  reflectance  to  space  by  clear  skies  (after 
Coulson,  1959)  could  be  effectively  approximated  by  least  squares  in 
the  following  form 

a (R)  =  .085  +  .25074  [log  (^-)  Sec  z) ]  (4-6) 

o 

where  P  =  1013.25  mb.  In  Eq.  (4-6),  tt/P  =  1  in  view  of  the  fact  that 
o  o 

mean  sea  level  pressure  7T  is  close  to  1000  mb.  Also 

Sec  z  =  (Cos  z) 

with  Cos  z  given  by  Eq.  (4-2). 

The  surface  albedo  a(G)  is  another  reflective  parameter  utilized  in 
this  study.  Over  oceanic  areas  the  following  formula  for  a(G)  after 
Gates  et  al  (1971) ,  was  utilized: 

a(G)  =  max  {.06,  .06  +  .54  (.7  -  Cos  z) }  .  (4-7) 

As  described  in  Section  III,  four  combinations  of  reference-cloud 
cases  are  possible  with  a  two- layer  cloud  model.  The  disposition  of 
F(S)  under  each  of  these  cases  will  be  discussed  in  the  remainder  of 
this  subsection. 

1.  Clear  Sky  Case 

In  the  clear  sky  (0,0)  case  the  F(S)  insolation  was  subjected 
to  both  Rayleigh  scattering  reflectance  a(R)  and  the  surface  reflectance 
a(G) .  Considering  the  likelihood  of  a  succession  of  multiple  reflections 
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between  earth  and  atmosphere,  the  F(S)  insolation  actually  penetrating 


the  earth’s  surface  after  scattering  is  given  by 


IS10(0,0)  =  F(S)  [l-a(R)]  [l+a(R)a(G)  +  ...  (a(R)a(G) ) 


n 


+  . ..]  * (l-a(G) ) 


(4-8a) 


that  is,  by 


isio(o,0)  =  F (s )  [l-a(R)  ]  [i-a(G)]/[i-a(R)a(G)  ] 


(4-8b) 


2 .  Cloudy-Sky  Cases 

In  the  three  cases  in  which  clouds  were  present,  F(S)  insolation 
absorbed  by  the  ground  at  each  gridpoint  was  computed  using  the  follow¬ 
ing  equation  (after  Arakawa,  1972) : 


IS10  (1 , 1)  =  F  (S )  (l-R(l))  (1-R(2)  )  (l-a(G)) 


*{l-[R(l)R(2)  +  R(2)a(G)  +  R(l)a  (G) 
-  2R(l)R(2)a(G)]  }_1  . 


(4-9) 


As  indicated  by  the  notation  (1,1),  denoted  CL^  =  CL^  =  1.0,  Eq.  (4-9) 
is  the  formula  used  in  calculating  F(S)  insolation  absorbed  by  the 
ground  in  the  case  where  overcast  clouds  are  present  at  both  levels  of 
Fig.  2.  Also  in  Eq.  (4-9),  initial  values  of  cloud-reflectance  were 
chosen,  namely  R(l)  =  .54  for  the  mid-level  clouds  between  k=4  and  6, 
and  R (2 )  =  .66  for  the  low-level  clouds  between  k=8  and  9.  Both  cloud- 
reflectance  values  are  as  suggested  by  C.  D.  Rodgers  (1967) . 

For  all  the  other  cloud  cases,  the  following  changes  were 
applied  to  Eq.  (4-9).  In  the  (1,0)  case  (CL^  =  1.0,  CL2  =  0.0),  the 
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desired  earth-absorbed  insolation  is  obtained  by  setting  R(2)  =0.0 
in  (4-9) ,  from  which  it  follows  that 

IS10 ( 1 , 0)  =  F(S)  (l-R(l)  )  (l-a(G)  )/[l-R(l)a(G)  ]  .  (4-10) 

In  the  case  (0,1),  one  sets  R(l)  =  0.0  in  (4-9)  so  that  (4-9)  simpli¬ 
fies  to 

IS10(0,1)  =  F(S)  (1-R (2 )  )  (l-a(G)  )/[l-R(2)a(G)  ]  .  (4-11) 

Note  that  with  a  cloud  overcast  present,  the  Rayleigh  clear-sky  scatter¬ 
ing  a(R)  does  not  appear  in  Eqs.  (4-9),  (4-10)  or  (4-11),  but  is  included 

empirically  in  the  cloud  reflectances  R(l)  and  R(2)  . 

3 .  Composite  F(S)  Insolation  at  Earth 

Equations  (4-8) ,  (4-9) ,  (4-10)  and  (4-11)  were  utilized  in  the 

computation  of  the  cloud-weighted  F(S)  insolation  penetrating  the 
earth 1 s  surface  considering  the  areal-weights  of  the  cloud  combinations 
denoted  by  (0,0),  (1,1),  (1,0)  and  (0,1)  about  a  gridpoint.  The  re¬ 

sultant  F (S )  insolation  penetrating  the  earth’s  surface  denoted  by 
IS10  is  therefore  expressible  as 

IS10(CL1,  CL2)  =  IS10(0,0)  W(0,0) 

+  IS10(1,1)  W(l,l) 

+  IS10(1,0)  W(1,0) 

+IS10(0,1)  W(0,1)  .  (4-12) 

Here  the  weighting  factors  W(0,0),  W(l,l),  W(1,0)  and  W(0,1)  are  com¬ 
puted  in  Eqs.  (2-3a,b,c,d)  respectively.  Note  finally  that  the  part 
of  F (S )  insolation  reflected  to  space  is  found  by  subtracting  IS10 
(CI^,  CL2)  from  F(S)  . 
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Table  VI  lists  the  F (S ) -disposition  resulting  from  a  particular  radia¬ 
tive  sounding  at  gridpoint  (1,1)  on  16  April  1974.  The  individual 
computations  of  IS10  as  they  apply  for  the  possible  overcast-clear 
layer  cases  are  made  under  the  heading  "ISIO."  The  difference 

REFS  =  F(S)  -  ISIO  (4-13) 

in  each  case  represents  F (S) -insolation  reflected  to  space  while 

STEAN  ■  -iSk  (4-14) 

has  been  computed  as  that  portion  of  the  F (S) -insolation  incident  at 
the  sea  surface  just  prior  to  transmission  by  the  surface.  Note  that 
no  absorption  in  air  has  been  included  in  the  computations  of  Table  VI, 
and  that  the  only  absorption  permitted  is  that  implicit  in  ISIO. 

Finally  at  the  bottom  of  each  column,  e.g.,  ISIO,  the  composite  value 
has  been  computed  by  means  of  the  weighting  scheme  of  Eq.  (4-12) . 

TABLE  VI.  A  sample  listing  of  values  of  F(S)  insolation  (ly  min 
computed  at  gridpoint  (1,1)  for  16  April  1974  using  Eqs.  (4-6,..., 

4-14) . 

Cloud-case  Weight 


(Cl  ,  cl2) 

W(CLlf  CL2) 

ISIO 

REFS 

STRAN 

(0,0) 

.0993 

.4305 

d751 

.5216 

(1,0) 

.1484 

.2539 

.3517 

.3076 

(0,1) 

.3016 

.1921 

.4135 

.2327 

(1,1) 

.4507 

.1340 

.4716 

.1696 

Composite-F (S) 

values 

.2014 

.4041 

.2441 
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C.  DISPOSITION  OF  F (A)  INSOLATION 


The  fractional  portion  of  the  solar  insolation  subject  to  absorption 
by  atmospheric  water-vapor  and  carbon  dioxide  is  covered  in  the  follow¬ 
ing  subsections. 

1.  Clear-sky  Case  (0,0) 

The  Manabe-Moller  absorptivity  function  provided  the  necessary 
absorptivity  values  for  the  key  layers  in  this  case.  The  form  of  this 
absorptivity  function  is 

303 

a(2,k)  =  .  271 [U (2 , k)  Sec  z] *  .  (4-15) 

Here  a  is  the  absorptivity  applied  to  the  pressure-scaled  water  vapor 
mass  between  levels  2  and  k  (Fig.  2)  along  the  zenith  slant-path  angle 
z.  The  resultant  absorbed  insolational  energy  in  the  particular  layer 
(2,4)  is  then  given  by  the  Manabe-Moller  relation 

303 

A24  =  0.271F(A)  [U(2,4)  Sec  z]  #  .  (4-16) 

In  the  same  manner  A26,  A28  and  A210  are  found.  Then  the  absorbed  in¬ 
solation  in  the  layers  (4,6),  (6,8)  and  (8,10)  are  computed  by 


A46  =  A26  -  A24 

(4-17) 

A68  =  A28  -  A26 

(4-18) 

A810=  A210-  A28  . 

(4-19) 

Water-vapor  mass  above  level  2  was  assumed  negligible  in  the  F (A)  dis¬ 
position  of  the  solar  insolation. 

By  subtracting  A210  from  F(A),  the  direct  transmission  of  F (A) 
insolation  impinging  at  the  earth's  surface  was  determined.  The  trans¬ 
mission  of  F (A)  insolation  is  then  further  reduced  by  the  transmissivity 
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(l-a(G))  after  surface-ref lectance, which  leads  to  the  earth-absorbed 
insolation 

IA10 (0,0)  =  F(A) {1-.271 [U(2, 10)Sec  z] *3°3} (l-a(G) ) .  (4-20) 

The  transmitted  energy  impinging  upon  the  earth  just  prior  to  absorp¬ 
tion  is 


TRANA  (0,0)  =  IA10(0,0)/[l-a(G)  ]  .  (4-21) 

In  the  remainder  of  this  subsection,  which  discusses  the  cloudy 
layer  cases,  representative  cloud  reflectivities  and  cloud  absorptivi- 
ties  were  initially  adopted,  after  C.  D.  Rodgers  (1967) ,  for  the  two 
possible  cloud  layers.  These  initial  cloud  reflectivities  were  RA(1)  = 
.46  and  RA(2)  =  .50  while  the  cloud  absorptivities  were  A(l)  =  .20  and 
A (2)  =  .30.  Note  that  the  reflectivities  for  the  F (A)  wavelengths  are 
somewhat  smaller  than  those  adopted  for  the  F(S)  wavelengths.  The  pro¬ 
cedure  of  considering  the  cloud  conditions  to  be  overcast  whenever  they 
appear  and  then  applying  the  appropriate  weighting  factors  in  the  com¬ 
posite  summation  will  again  be  followed  as  in  Sec.  IV. B. 

To  simplify  the  discussion  for  the  cloud-covered  cases,  the  (1,1) 
case  will  be  presented  first  as  it  contains  representative  type  equa¬ 
tions  for  the  remaining  two  cases,  (1,0)  and  (0,1). 

2 .  Overcast  in  Both  High-  and  Low-cloud  Layers 

The  following  set  of  formulas  illustrate  the  model  disposition 
of  incoming  solar  insolation  (F(A))  from  level  k=2  to  the  earth* s  sur¬ 
face  and  permits  determination  of  the  amount  of  insolation  absorbed  by 
the  atmospheric  layers  and  by  the  earth* s  surface.  The  dashed  separa¬ 
tion  lines  are  introduced  to  subdivide  the  absorption  and  reflection 
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physics  of  the  model  into  subsections  which  permit  the  analysis  to 
proceed  more  or  less  within  successive  200  mb  layers 0  The  equations 
relate  to  the  parameters  in  Fig*  6  where  the  insolations ,  A24 ,  A46, 

A68,  A89  and  A910,  etc.,  represent  the  contributions  to  the  insola¬ 
tion  absorbed  in  the  layers  involved.  Symbols  F2,  F4,  F6,  F8  and  F9, 
etc.,  depict  the  streams  of  insolation  passing  through  the  indicated 
level.  A  vertical  arrow  implies  the  direction  of  insolation  passage, 
i.e.,  1  denotes  downward  insolation,  f  upward-reflected  insolation, 
and  downward-reflected  insolation.  Terms  involving  the  symbols 
"TD",  as  expressed  by  the  functions  of  Eq.  (4-22),  indicate  the  Manabe- 
Moller  transmissivities  for  diffuse  insolation  in  the  layer  beneath  an 
existing  cloud.  In  the  latter  situation,  the  term  Sec  Z  in  Eq.  (4-15) 
is  effectively  replaced  by  the  mean  slant-path,  Sec  Z  —  5/3  (Katayama, 
1966) ,  e .g . , 

TD68  =  1- . 271 [U (6,8)  5/3] *3°3 

303 

TD910  =  1-.271 [U (9,10) 5/3] *  .  (4-22) 

In  accordance  with  the  above  definitions  and  as  described  by  Fig.  6, 
the  formulas  for  the  (1,1)  case  are  listed: 

F4+  =  F(A)  (1-. 271  [U (2,4)  Sec  z]*303) 

F4+  =  F44'  (RA  (1) ) 

303 

F2+  =  F4+  (1- . 271 [U (2 , 4)  Sec  z]  ) 

A24  =  F  (A)  -  F44-  +  F4+  -  F2+  (4-23a) 
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A46  =  F4+(A(1)) 


(4-23b) 


F64  = 

F41  -  F4t  -  A46 

F84  = 

F64 (TD68) 

F8t  = 

F84  (RA (2) ) 

F6t  = 

F8t (TD68) 

F644= 

F6t (RA (1) ) 

F844= 

F644  (TD68) 

A68  = 

F6l  -  F84  +  F8t  -  F6t  +  Fell  -  F844  (4-23c) 

A89  = 

(F84  +  F844)  (A ( 2)  ) 

F94-  = 

F84  -  F8t  +  F844  -  A89 

F104= 

F94 (TD910) 

F10t= 

FI 04  (a(G) ) 

F9t  = 

FlOt (TD910) 

F9i4= 

F9t (RA (2) ) 

F10H=  F9H  (TD910) 

A910  =  F9I  -  F104  +  FlOt  -  F9+  +  F944  -  F1044 

A810  =  A89  +  A91C  (4-23d) 

Note  that  the  effect  of  multiple  reflections  between  clouds  or 
between  the  earth's  surface  and  the  lower  cloud  has  been  incorporated 
to  include  only  the  effect  of  two  reflections,  with  the  lowermost  re¬ 
flecting  surface  absorbing  the  remaining  impinging  insolation.  Computa 
tions  indicated  that  the  insolation  remaining  after  two  reflections  was 
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Figure  6.  Schematic  representation  of  F(A)  insolation  disposition  in 
the  case  of  two  overcast  layers. 


50 


too  small  to  warrant  the  consideration  of  further  reflections.  Also 
insolation  reflected  upward  from  a  lower  interface  (cloud  or  ground) 
to  the  base  of  an  upper  cloud  deck  has  not  been  subjected  to  upper 
cloud  absorption.  This  tends  to  reduce  very  slightly  the  secondary 
cloud-absorption. 

From  Eq.  (4-23d) ,  the  impinging  F (A)  insolation  at  the  earth’s 
surface  can  be  expressed 


TRANA (1,1)  =  FlOl  +  F10+  + 


(4-24) 


The  F (A)  insolation  which  is  actually  absorbed  by  the  earth's  surface 
(see  Fig.  6)  may  be  written  as 


IA10 (1,1)  =  F10l(l-a(G))  +  F10++ 


(4-25) 


3 .  Disposition  of  F(A)  Insolation  with  an  Upper  Overcast  Only 

With  a  single  cloud  layer  present  only  at  the  upper  level  (Fig. 
7)  the  equations  depicting  the  model-disposition  of  incoming  insolation 
becomes  a  simplified  subset  of  the  previous  case: 


F4+  =  F (A)  (1- . 271  [U ( 2 , 4)  Sec  z] 


.303 


F4t  =  F4+ (RA ( 1) ) 

F2+  =  F4t(l-.271[U<2,4)  Sec  z]*303 


A24  =  F (A)  -  F4+ +  F4i  -  F2t 


(4-26a) 


A46  =  F44-  (Al)  ) 


(4-26b) 
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Figure  7. 


Schematic  representation  at  F(A)  insolation  disposition 
with  an  upper  overcast  layer  only. 
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F64- 

F44-  -  F4t  -  A46 

F84- 

F64-  (TD68) 

F104-  = 

F64-  (TD610) 

F  lOt  = 

F104-  (a(G)  ) 

F8+ 

FlOt (TD810) 

F6f 

F10I (TD610) 

F64-+  = 

F6t (RA ( 1) ) 

F84-4-  = 

F64-4-  (TD68) 

FlOii  = 

F84-4-  (TD810) 

A68 

F64-  -  F8+-  +  F8t  -  F6+  +  F6ll  -  F844  (4-26c) 

A810  = 

F84-  -  F10+  +  FlOf  -  F8+  +  F844  -  F1044  (4-26d) 

The  variables  used  above  are  defined  in  a  similar  manner  to  those  in  the 
(1,1)  case*  The  impinging  insolation  at  the  earth’s  surface  may  be  for¬ 
mulated  as  follows 

TRANA  (1,0)  =  F10+  +  FlOK  (4-27) 

while  the  F(A)  insolation  absorbed  by  the  earth  in  this  case  is  given  by 

IA10  (1,0)  =  Fl0i(l-a(G))  +  F1044  .  (4-28) 
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4.  Disposition  of  F (A)  Insolation  with  a  Low  Overcast  Only 


With  an  overcast  lower  cloud  layer  the  model-disposition  symbols 
are  as  depicted  in  Fig.  8  and  are  physically  related  as  follows: 


F44- 

=  F  (A)  (1-  .  271  [U  (2 ,4) 

_  ,  •  303 

Sec  z]  ) 

F64- 

=  F  (A)  (1- .  271  [U  (2 , 6) 

_  ,  -  303 

Sec  z]  ) 

F84- 

=  F (A) (1- . 271 [U ( 2 , 8) 

o  , .303, 

Sec  z]  ) 

F8t 

=  F8HRA(2)  ) 

F6t 

=  F8t(l-.271 [U  (6 , 8) 

,  .303, 

Sec  z]  ) 

F4  + 

=  F8+  (1-  .  271  [U  (4 ,8) 

„  ,-303, 

Sec  z]  ) 

F2  + 

=  F8t (1- . 271 [U(2,8) 

„  ,.303, 

Sec  z]  ) 

A24 

=  F  (A)  -  F44-  +  F4  +  - 

F2  + 

(4-29a) 

A46 

=  F44-  -  F6+  +  F6+  - 

F4  + 

(4-29b) 

A68 

=  F64-  -  F84-  +  F8I  - 

F6t 

(4-29c) 

A89 

=  F84-  (A  ( 2)  ) 

F94- 

=  F84-  -  F8+  -  A89 

F104- 

=  F94-  (TD910) 

F10+ 

=  Fl(H  (a  (G) ) 

F9i 

=  fio-Mtosio) 

F94-4- 

=  F9+ (RA(2) ) 

FlO'1'4- 

=  F94-4-  (TD910) 

A910 

=  F94-  -  F104-  +  FlOt 

-  F9+  +  F94-4-  -  F1044 

A810 

=  A89  +  A910 

(4-29d) 
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Figure  8.  Schematic  representation  of  F(A)  insolation  disposition 
with  a  lower  overcast  layer  only. 
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Again  the  variables  are  defined  as  stated  before  in  the  (1,1) 


case.  Likewise  the  incident  flux  is  defined  at  the  earth* s  surface  as 


TRANA (0,1)  =  FlOl  +  FlOll 


(4-30) 


while  that  portion  which  is  absorbed  by  the  earth  (Fig.  8)  is 


IA10 (0,1)  =  Fl0i(l-a(G))  +  F1011 


(4-31) 


In  Eqs.  (4-25) ,  (4-28) ,  and  (4-31) ,  the  quantity  FlOil  is  small  enough 
in  each  case  so  that  no  further  reflections  from  the  earth  were  con¬ 
sidered. 

5.  Composite  F(A)  Layer-Absorptions  and  Surface-Absorption 

Insolation 

As  has  been  previously  discussed,  the  standard  grid-area  weight¬ 
ing  scheme  of  this  study  was  applied  to  obtain  composite  values  of  the 
absorbed  F (A) -insolation  in  key  layers  and  also  within  the  earth  * s  sur¬ 
face.  The  weighting  factors  applied  to  the  corresponding  overcast- 
combination  absorption  quantities  provided  the  following  composite 
results : 


A24(CL1,CL2)  =  A24  (0,0)  W(0,0)  +  A24(l,0)  W(1,0) 


+  A24 (0,1)  W(0,1)  +  A24 (1,1)  W(l,l) 


(4-32) 


A46(CL1,CL2)  =  A46  (0,0)  W(0,0)  +  A46(l,0)  W(1,0) 


+  A46  (0,1)  W(0,1)  +  A46  (1,1)  W(l,l) 


(4-33) 


A46(CLl,CL2)  =  A68 (0 , 0)  W(0,0)  +  A68(l,0)  W(1,0) 


+  A68  (0,1)  W(0,1)  +  A68  (1,1)  W(l,l) 


(4-34) 
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(4-35 ) 


A810(CL1/CL2)  =  A810 (0, 0)  W(0,0)  +  A810(l,0)  W(1,0) 

+  A810  (0,1)  W(0,1)  +  A810)  1,1)  W(l,l) 

IA10  (CL1  ,  CL0 )  =  IA10  (0,0)  W(0,0)  +  IA10(1,0)  W(1,0) 

+  IA10(0,1)  W  (o ,  1)  +  IA10  (1,1)  W(l,l)  •  (4-36) 

The  weighting  factors  W (0,0) , . • . , (Wl, 1)  were  listed  in  Eqs .  (2-3a ,b, c ,d) , 

and  A24 (0,0) ,  A46(0/0)/  A68(0,0),  A810(0,0)  and  IA10(0,0)  are  given  in 
each  clear  sky  case  (0,0)  about  each  gridpoint  by  Eqs.  (4-15)  ,  (4-16)  , (4-17) 
(4-18) ,  (4-19)  and  (4-20)  respectively. 

The  results  for  the  absorption  in  layers  (2,4) ,  (4,6) ,  (6,8) , 

(8,10)  and  at  the  surface,  level  k  =  10,  are  shown  in  the  following  table: 

TABLE  VII.  A  sample  listed  of  F (A)  insolation  values  (ly  min  1) 


computed 

at  gridpoint 

(1,1)  for 

the  16 

April 

case. 

Cloud-Case  Weight 

CL  ,  CL2  W(CLlfCL2) 

A24 

A46 

A68 

A810 

IA10 

REFA 

TRANA 

(0,0) 

.0993 

.0454 

.0356 

.0386 

.0418 

.1348 

.0285 

.1634 

(1,0) 

.1484 

.0633 

.0559 

.0280 

.0156 

.0482 

.1137 

.0581 

(0,1) 

.3016 

.0464 

.0381 

.0728 

.0771 

.0239 

.0664 

.0286 

(1,1) 

.4507 

.0633 

.0559 

.0432 

.0285 

.0104 

.1233 

.0126 

Composite 

-F (A)  values 

.0564 

.0485 

.0494 

.0425 

.0324 

.0953 

.0391 

In  the  computational  scheme  indicated  by  the  entries  of  Table 
VII,  the  reflected  F (A)  insolation  to  space  has  been  depicted  by  the 
symbol  REFA,  and  its  values  follow  from 

REFA  =  F (A)  -  A24  -  A46  -  A68  -  A810  -  IA10  (4-37) 
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whereas  the  TRANA  dispositions  are  given  by  Eqs .  (4-21),  (4-24),  (4-27) 


and  (4-30)  or  by  its  weighted-mean  value  in  the  case  of  TRANA-composite . 

6,  Absorptivity  (ABA)  by  Layers 

Here  the  (fractional)  absorptivity  as  well  as  the  actual  insola¬ 
tion  values  absorbed  in  the  layers  are  considered.  In  the  computation 
of  absorptivity,  which  is  fractional  absorption,  the  total  undepleted 
insolation  at  the  top  (k=0)  is  used  as  a  base.  The  following  equation 
was  utilized  in  this  calculation: 

-2 

FADJ  =  2.00  (r/r  )  Cos  z  .  (4-38) 

m 

The  absorptivity  of  the  troposphere  ABA  was  computed  from  the  ratio  of 
the  insolation  absorbed  in  the  troposphere  to  the  insolation  incident 
at  the  top  of  the  atmosphere  rather  than  at  k=2 : 


ABA 


A24  +  A46  +  A68  +  A810 
FADJ 


(4-39) 


D.  ALBEDO  (ALB)  OF  THE  EARTH -TROPOSPHERE  SYSTEM 

In  considering  the  planetary  albedo,  the  reflected  insolations  of 
the  earth-troposphere  systems  in  both  the  F(A)  and  F(S)  insolational 
regions  must  be  recalled  by  the  program.  Thus  REF  is  computed  at  each 
gridpoint  as  the  sum  of  the  reflected  insolation  energy  in  F(A),  denoted 
REFA  in  Eq.  (4-37)  and  the  reflected  part  of  F(S)  denoted  REFS  in  Eq. 
(4-13) : 


REF  =  REFS  +  REFA  (4-40) 

Finally  the  planetary  albedo  is  related  to  FADJ  through 
REF 


ALB 


'FADJ  # 


(4-41) 


E.  COMPOSITE  ABSORPTIVITY  (ABG)  BY  THE  EARTH- SURFACE ; 

COMPOSITE  ATMOSPHERIC  TRANSMISSIVITY  (ATRAN) 

1.  Absorptivity  (ABG)  of  Earth 

By  summing  the  weighted  values  of  F(S)  and  F(A)  portions  of 
the  incoming  insolation  entering  the  earth,  the  total  insolation  ab¬ 
sorbed  at  the  earth's  surface  was  computed.  This  quantity  when 
divided  by  the  extraterrestrial  insolation  gave  the  fractional  absorp¬ 
tivity  (ABG)  of  the  earth's  surface.  The  equation  for  ABG  was 


ABG  = 


IA10  +  IS10 
FADJ 


(4-42) 


where  IA10,  IS10,  and  FADJ  were  defined  previously  by  Eqs.  (4-36), 
(4-12)  and  (4-38)  respectively. 

2 .  Transmissivity  (ATRAN)  of  the  Troposphere 

Also  computed  was  the  total  insolational  energy  TRAN,  incident 
at  the  earth's  surface  just  before  absorption  by  the  surface.  This 
calculation  is  given  by 


TRAN  =  TRANA  +  STRAN  . 


(4-43) 


Here  STRAN  =  [IS10/ (1-a (G) ]  was  previously  defined  in  Eq.  (4-14)  and 
a(G)  was  given  in  Eq.  (4-7) .  TRANA  has  also  been  defined  as  the  weighted 
value  of  TRANA  (0,0),  TRANA  (1,1),  TRANA  (1,0)  and  TRANA  (0,1)  given  by 
Eqs.  (4-21),  (4-24),  (4-27)  and  (4-30).  Note  also  in  justification  of 

STRAN  that  the  four  cases  for  IS10  of  (4-8),  (4-9),  (4-10),  and  (4-11) 

each  have  the  common  factor  (l-a(G))  in  the  numerator  and  therefore  each 
transmitted  F(S)  insolation  component  available  at  the  earth  just  before 
absorption  needs  only  be  divided  by  (l-a(G)).  TRAN  may  thus  be  viewed 
as  the  total  insolational  energy  incident  at  a  pyrheliometer  located  at 
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earth.  The  (fractional)  transmissivity  of  the  troposphere  (ATRAN)  is 


then  computed  from 


ATRAN  E  TRAN/FAD J  .  (4-44) 

Note  finally  that  the  major  dispositions  of  the  total  insolation  at 
the  indicated  map  times  have  now  been  identified  by  the  fractional 
values,  ALB,  ABA  or  ABG,  and  ATRAN,  representing  the  reflectivity 
(albedo) ,  absorptivity  of  air  or  earth,  or  atmospheric  transmissivity 
as  the  case  may  be. 

F.  ALBEDO  TUNING  BY  COMPARISONS  WITH  SATELLITE  CLIMATOLOGY 
1 .  General  Remarks  Concerning  a  Need  for  Tuning  Albedos 

In  the  four  previous  radiational  studies  (A,  B,  C,  D)  compari¬ 
sons  were  made  between  the  solar-insolation  albedo-model  (ALB MOD)  com¬ 
putations  and  the  satellite-climatology  albedo  (ALBRAS)  of  Raschke 
et  al.  (1973).  The  results  of  these  comparisons  indicated  excessively 
high  values  of  ALBMOD  especially  in  the  tropical  and  subtropical  oceanic 
areas.  The  question  was  raised  by  the  previous  investigators  as  to 
whether  vertically-structured  convective  cloud  elements  in  the  tropics 
would  have  as  high  a  reflective  capability  as  attributed  to  the  large- 
scale  cloud  masses  existing  primarily  in  horizontal  layers  as  indicated 
by  the  parameterization  of  Eq.  (2-2) .  It  should  be  recalled  that  the 
initial  cloud-reflectances  used  were  as  suggested  by  C.  D.  Rodgers 
(1967)  and  are  listed  in  Secs.  IV. B. 2.  and  IV.C.l. 

For  the  purpose  of  tuning  the  model  solar  dispositions,  and 
in  particular  ALBMOD  with  ALBRAS,  the  initial  values  of  Rodgers*  cloud- 
reflectances  were  simply  adjusted  by  a  multiplicative  factor  (f)  which 
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turns  out  by  the  least-squares  fitting  to  be  smaller  than  unity.  This 
device  allows  a  greater  fraction  of  solar  radiation  to  penetrate  down¬ 
ward  through  the  cloud-types  in  both  the  tropical  and  subtropical  areas 
and  provides  better  agreement  between  ALBMOD  and  ALBRAS. 

2 .  Method  of  Tuning  ALBMOD  to  ALBRAS 

The  model- tuning  process  involved  several  steps.  First,  the 
new  variable  f  was  defined  and  allowed  to  range  in  value  from  1.2  to 
0.25  by  0.05  intervals.  The  f-values  were  then  multiplied  by  the 
initially  used  cloud-reflectance  values  (after  Rodgers,  1967)  giving 
twenty  sets  of  new  cloud-reflectivity  values.  For  each  set  of  new 
values,  ALBMOD  and  subsequently  the  difference  A (ALB)  between  ALBMOD 
and  ALBRAS ;  i . e . , 


A (ALB)  =  ALBMOD  -  ALBRAS  (4-45) 

was  recomputed  at  each  gridpoint  for  the  four  midseasonal  dates. 

Because  cloud-types  differ  in  structure  over  tropical  and  extra 
tropical  oceanic  areas  the  resulting  A(ALB)  data,  from  (4-45),  were 
divided  into  two  populations  denoted  TR  and  HL.  The  population  TR  in¬ 
cluded  the  set  of  gridpoints  located  south  of  and  including  25N,  while 
HL  included  the  set  of  gridpoints  located  north  of  25N  over  all  four 
meridians . 

Only  those  gridpoints  where  the  value  of  at  least  one  of  CL^ 
and  CL^  was  greater  than  0.1  were  included  in  those  populations.  The 
gridpoints  so  excluded  were  considered  climatologically  unrepresenta¬ 
tive  for  the  purpose  of  cloud- tuning  the  albedo.  The  A (ALB)  data  was 
additionally  divided  by  season  giving  a  total  of  eight  separate  cases 
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further  defined  in  sequence  as  WIN-TR,  WIN-HL,  SPR-TR,  SPR-HL,  SUM-TR, 


SUM-HL,  AUT-TR  and  AUT-HL. 

The  Bimedical  set  of  programs  (Dixon,  1973)  was  utilized  to 
compute  the  root-mean-square  (RMS)  of  A (ALB)  values  corresponding  to 
a  wide  range  of  f-values  for  each  of  the  seasons  and  regions.  Values 
of  f  leading  to  the  minimum  RMS  of  A (ALB)  were  isolated  and  appear  in 
Table  VIII,  along  with  the  RMS  difference  of  A (ALB). 


TABLE  VIII.  Values  of  f  and  corresponding  minimum  RMS 
values  of  A (ALB). 


f 

Minimum  RMS 
of  A (ALB) 

WIN-TR 

.40 

.069 

WIN-HL 

• 

00 

o 

.073 

SPR-TR 

.35 

.067 

SPR-HL 

.65 

.082 

SUM-TR 

.60 

.079 

SUM-HL 

.70 

.109 

AUT-TR 

.40 

.052 

AUT-HL 

.60 

.073 

As  can  be  seen  from  Table  VIII,  the  listed  values  of  f  in 
winter  indicates  the  largest  difference  between  the  tropical  and  extra- 
tropical  cloud- types.  The  evenly  layered  cloud  formations  of  the 
northern  latitudes  in  mid-January  seems  more  closely  to  approximate 
those  considered  by  C.  D.  Rodgers  (1967).  On  the  other  hand,  the 
smallest  difference  in  f-values  occurs  between  SUM-TR  and  SUM-HL 
indicating  a  less  significant  difference  in  cloud  types  over  the 
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complete  latitude  range  during  mid-July.  As  would  be  expected  the 
Spring  and  Autumn  f-ranges  fall  in  between  the  Winter  and  Summer  cases. 

A  full  listing  of  zonally-averaged  ALBMOD  values,  after  tuning  at  all 
gridpoints,  and  of  zonally-averaged  ALBRAS  values,  generated  from 
identical  gridpoints,  can  be  found  in  Table  IX. 

Upon  examination  of  Table  IX  it  is  apparent  that  the  tuned 
model-albedo  values  compare  closely  with  the  satellite  albedos  derived 
from  Raschke.  Table  IX  includes  the  annual  averages  and  the  Northern 
Hemisphere  cosine-weighted  means  as  well  as  the  mid-seasonal  averages. 

Further  inspection  of  the  weighted  means  shows  that  ALBMOD  is 
somewhat  smaller  than  ALBRAS  for  each  of  the  four  seasons  and  annual 
average.  However  most  of  this  underestimate  seems  to  occur  in  the 
HL-population  where  the  result  of  minimizing  A (ALB)  led  to  values  of  f 
substantially  less  than  unity  (Table  VIII) .  This  result  was  unexpected 
and  is  probably  due  to  a  mismatch  between  the  model  clouds  (Eq.  (2-2)) 
and  those  averaged  in  the  satellite  climatology  particularly  in  the 
higher  latitudes. 

The  conclusion  is  that  reflectance-tuning  to  satellite  albedo 
appears  valid  in  the  tropics  and  subtropics  where  cloud-formations 
on  a  given  day  tend  to  be  persistent.  However  tuning  by  this  method  in 
high  latitudes  is  less  conclusive  because  interdiurnal  cloud  variability 
is  much  greater. 

In  any  case,  all  further  computations  of  solar  dispositions  in 
this  study  will  utilize  the  appropriate  cloud-reflectance  tuning  factors 
as  listed  in  Table  VIII,  depending  upon  season  and  latitude. 
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V.  MERIDIONAL  CROSS-SECTIONAL  DEPICTION 
OF  THE  RADIATIVE-BALANCE  COMPUTATIONS 


A.  GENERAL 

The  general  design  of  this  section  is  to  utilize  all  of  the  compu¬ 
tational  concepts  discussed  in  Sections  III  and  IV  in  the  computations 
for  a  single  time-step  in  the  radiative  heating  package  developed  for 
use  in  the  FNWC  prediction  model.  After  testing  ALBMOD  with  the  corres¬ 
ponding  values  of  Raschke  et  al.  (1973) ,  it  was  decided  that  only  the 
computations  by  the  appropriate  cloud-reflectance  tuning  factors  (Table 
VIII)  would  be  displayed  in  the  meridional  cross-sections  (Figs.  10... 
17)  and  in  the  comparisons  of  the  four-layer  versus  the  two-layer  cool¬ 
ing  rates  summarized  in  Table  X. 

The  radiative  calculations  were  performed  at  each  gridpoint  of  the 
four  meridians  for  each  of  the  four  mid-seasonal  days  considered;  how¬ 
ever,  only  the  winter  and  summer  cases,  as  illustrating  the  two  extreme 
situations,  are  presented  in  this  section. 

B.  GEOGRAPHICAL  REPRESENTATION  OF  THE  RADIATIVE -BALANCE  DISTRIBUTION 

The  FNWC  gridpoint  processed  analyses  for  0000GMT  were  used  at  the 

three  Pacific  cross-sections,  while  that  for  1200GMT  were  used  for  the 
single  Atlantic  meridian.  This  was  done  so  that  the  set  of  gridpoints 
was  considered  to  be  subject  to  the  actual  radiative- transfer  calcula¬ 
tions  involved  for  these  specific  times.  Figure  9  depicts  in  symbolic 
language  the  key  to  the  computational  entries  in  Figs.  10,  11,... 16, 

17.  This  symbolic  list  presents  the  computations  made  at  each  radia¬ 
tive  sounding  gridpoint  ( I, J) having  data  in  the  form  of  Table  I.  The 
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computations  proceed  from  the  tropopause  (approximately  level  k=2)  to 
the  ocean  surface.  For  purposes  of  climatological  data  comparison, 

Figs,  10,  11,. ..16,  17  were  developed  by  interpolating  gridpoint  results 
to  integral  multiples  of  5-degree  latitudinal  increments.  The  inter¬ 
polation  routine  to  this  gridpoint  spacing  made  use  of  the  Lagrangian 
cubic  interpolation  scheme  (after  Spaeth,  1975) . 


Q(D  = 


(1-2)  (1-3)  (1-4) 
(1-2)  (1-3)  (1-4) 


O  (I-D  d-3)  (1-4) 
y2  (2-1)  (2-3)  (2-4) 


(1-1)  (1-2)  (1-4) 
y3  (3-1)  (3-2)  (3-4) 


(1-1) (1-2) (1-3) 
~4  (4-1)  (4-2)  (4-3) 


(5-1) 


Finally  for  ease  in  reconciling  the  magnitudes  of  all  radiative-trans¬ 
fer  rates,  the  time-dependent  solar  disposition  rates  have  been 
averaged  to  24-hourly  rates . 


C.  EXPLANATION  OF  SYMBOLIC  TERMS 
1 .  Cross-Section  at  Level  k=2 

The  discussion  of  all  insolation  parameters  discussed  previous¬ 
ly  in  Section  IV  dealt  with  the  specific  time  of  day  that  corresponded 
to  the  hour  angle  h  for  the  instantaneous  time  t  under  consideration. 
The  incident  solar  insolation  dealt  with  is  then 


F  (2)  =  S(— )  2  Cos  z  . 
r 

m 


(5-2) 


In  order  to  avoid  reference  to  specific  map  times  t,  the  instantaneous 
solar  hour-angles  were  h  =  35,  10,  55,  and  35  degrees,  respectively 
for  cross-sections  1,2,3  and  4  as  depicted  in  Figs.  10,  11,  12  and  13 
for  winter,  and  in  Figs.  14,  15,  16,  and  17  for  summer. 
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QAVE  represents  the  24-hour  average  of  F(2)  and  appears  as  the 
first  input  symbol  in  Fig.  9.  Its  value  is  considered  to  be  more  re¬ 
presentative  climatologically  for  the  data  day  under  consideration  than 
F  (2)  . 

QAVE  is  derived  by  the  formula 

«RVE  ■ F<2)  §sri  <5-3) 


where 

Cos  z  =  [H  Sin(f)Sin6  +  Cos4>Cos6SinH]  /TT  (5-4) 

H  =  ArcCos  [-Tan<{>Tan<5]  (5-5) 

Here  6  is  the  appropriate  solar  declination  angle  as  listed  in  Table  V, 
and  H  is  the  appropriate  hour  angle  at  local  sunset  at  latitude  (j).  The 
value  of  H  also  depends  upon  which  mid-seasonal  date  is  being  considered. 
Cos  z  in  Eq.  (5-4)  is  equal  to  the  24-hour  average  cosine  of  the  zenith 
angle,  Eq.  (4-2) .  The  24-hour  time  averaging  period  for  QAVE  gives  heat¬ 
ing  results  consistent  in  magnitude  with  the  terrestrial  flux  divergences, 
which  change  only  slightly  with  the  time  of  day.  The  conversion  to  ex¬ 
pected  daily  averaged  solar  disposition  quantities  is  compatible  with 
the  determination  of  a  hemispheric  radiative  balance  for  the  given  mid- 
seasonal  dates  (Sec.  IV. A.). 

Other  parameters  needed  for  level  k=2  are 

pr-)  c  7 

QREF  =  REF  (t)  ( - )  (5-6) 

Cos  2 

where 

REF (t)  =  F (2)  -  A24  -  A46  -  A68  -  A810 

-  (IA10  +  IS10)  (5-7) 


67 


REF ( t)  is  the  instantaneous  solar  reflected  insolation  at  a  gridpoint 
and  QREF  is  its  24-hour  average,  assuming  that  the  instantaneous 
planetary  albedo  remains  constant  for  the  24-hour  period.  This  assump¬ 
tion  requires  that  the  cloud  amounts  computed  at  the  indicated  synoptic 
times  are  representative  of  the  entire  day. 

The  same  principle  will  be  used  with  regard  to  all  other  solar 
parameters  in  the  conversion  from  time-dependent  values  at  solar  time  t 
to  24-hour  averaged  values.  Superior  bars  (  )  are  not  used  in  the  sym¬ 

bolism  for  the  averaged  values  shown  in  the  cross-sections  key.  Fig.  9, 
but  are  implied  by  the  use  of  QAVE,  etc.  in  the  solar-disposition  terms. 
The  24-hour  average  tropopause  balance,  BALT,  is  computed  from 

BALT  =  QAVE  -  (QREF  +  F2*)  (5-8) 

for  the  level  k=2  at  the  indicated  latitude.  Net  terrestrial  fluxes, 

* 

such  as  F^  ,  were  considered  to  be  constant  throughout  the  24-hour  period, 
a  valid  assumption  if  the  cloud  cover  remains  quasi-constant  for  the 
period. 

2 .  Cross-Sections  in  Layers  (2 , 4) ,  (4 , 6)  , (6 , 8)  and  (8,10) 

The  following  definitions  apply  for  the  four  layers  identified 
in  the  present  section  heading  and  are  further  identified  in  Fig.  9. 

All  of  the  heat  transfers  shown  in  these  layers  are  assumed  to  be  of 
radiative  character  only.  The  daily-averaged  radiative  heating  (cooling) 
rate  in  layer  (2,4)  is  given  by 

BAL  24  =  024  -  F24  (5-9) 
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Figure  10(a).  125W  Longitudinal  cross-section,  tropical  section.  Refer  to  Fig.  9  for  key.  Values 

computed  from  data  for  16  January  1974. 
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10(b)'.  12 5W  Longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key 

computed  from  data  for  16  January  1974. 
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11(b).  170W  longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key 

computed  from  data  for  16  January  1974. 
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Figure  12(a).  145E  Longitudinal  cross-section,  tropical  section.  Refer  to  Fig.  9  for  key.  Values 

computed  from  data  for  16  January  1974. 
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12(b).  145E  Longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key. 

computed  from  data  for  16  January  1974. 
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Figure  13(a).  35W  Longitudinal  cross-section,  tropical  section.  Refer  to  Fig.  9  for  key.  Values 

computed  from  data  for  16  January  1974. 
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13(b)0  35W  Longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key 

computed  from  data  for  16  January  1974. 
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Figure  14(a).  125W  Longitudinal  cross-section,  tropical  section.  Refer  to  Fig.  9  for  key.  Values 

computed  from  data  for  16  July  1974. 
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14(b).  125W  Longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key 

computed  from  data  for  16  July  1974. 
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15(b).  170W  longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key 

computed  from  data  for  16  July  1974. 
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Figure  16(a).  145E  Longitudinal  cross-section,  tropical  section.  Refer  to  Fig.  9  for  key.  Values 

computed  from  data  for  16  July  1974. 
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16(b).  145E  Longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  key 

computed  from  data  for  16  July  1974. 
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Figure  17(a).  35W  Longitudinal  cross-section,  tropical  section.  Refer  to  Fig.  9  for  key.  Values 

computed  from  data  for  16  July  1974. 
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17(b).  35W  Longitudinal  cross-section,  higher  latitude  section.  Refer  to  Fig.  9  for  kev 

computed  from  data  for  16  July  1974. 


where 


Q 24  =  daily-averaged  solar  absorption  in  layer  (2,4) 

and  is  defined  relative  to  A24(t)  by  a  cosine  transformation  similar 
to  Eq.  ( 5—3 ) ,  and 

F24  =  terrestrial  cooling  rate  in  (2,4)  . 

Similarly  the  daily-averaged  radiative  heating  (cooling)  rates  in 
layers  (4,6),  (6,8)  and  (810)  are  given  by 


BAL46  =  Q46  - 

F46 

(5-10) 

BAL68  =  Q68  - 

F68 

(5-11) 

BAL810=  Q810- 

F810 

(5-12) 

respectively . 

3 .  Cross-Section  at  Air-Sea  Interface  (k=10) 

The  radiative  balance  at  the  earth* s  surface  (BALB)  is  as 
defined  in  the  following  equation: 

•k 

BALB  =  QA3G  -  F  .  (5-13) 

QABG  is  the  24-hour  average  solar  insolation  absorbed  by  the  surface 
as  follows: 

QABG  =  QABG (t) (Cos  z/Cos  z)  .  (5-14) 

D.  MERIDIONAL  CROSS-SECTIONS  OF  THE  VERTICAL  RADIATION  BALANCE 
Figs.  10,  11,... 16,  17  as  previously  explained  represent  the 
single  time  step  of  heating  computations  for  each  of  the  four  meridians 
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used  in  this  study.  The  eight  figures  have  been  divided  into  (a)  tropi¬ 
cal  results  and  (b)  mid-to-high-latitude  results  for  the  mid-winter  and 
mid-summer  cases  respectively.  While  the  results  depicted  in  these 
cross-sections  are  exhibited  as  representing  daily-averaged  values,  they 
are  actually  based  upon  radiation-computations  at  the  specific  map  times 
of  OOOOGMT  and  1200GMT  on  16  January  and  16  July  1974.  Therefore,  for 
these  results  to  be  meaningful  as  a  stepwise  part  of  the  heat  package 
subroutine  of  FNWC,  the  solar  radiative  absorption  and  reflectance  terms 
would  have  to  be  recoverable  as  a  function  of  GMT,  i.e., 

F ( 2 , t)  =  QAVE*  (Cos  z/Cos  z)  (5-15) 

REF ( t)  =  QREF*  (Cos  z/Cos  z)  (5-16) 

etc.  Thus  solar  disposition  terms  may  then  be  utilized  in  connection 
with  the  one-hour  stepwise  application  of  the  thermodynamic  equation  of 
the  set  of  primitive  equations  used  in  the  FNWC  prediction  process, 
assuming  the  "2/3-CL"  parameterization  of  Eq.  (2-2) . 

E.  COMPARISON  OF  THE  FOUR- LAYER  WITH  THE  TWO-LAYER  FLUX  DIVERGENCE  MODEL 
In  previous  studies  (A,  B,  C,  D)  which  considered  only  the  two-layer 
flux-divergence  model  [layers  (2,6)  and  (6,10)]  it  was  assumed  that  the 
radiative-cooling  rate  in  layer  (2,6),  BAL26,  was  uniformly  distributed 
over  layers  (2,4)  and  (4,6).  A  similar  remark  is  applicable  to  BAL610 
relative  to  the  sublayers  (6,8)  and  (8,10).  In  an  effort  to  determine 
the  vertical  resolution  of  the  atmospheric  flux-divergences  separately, 
the  four  values  of  BAL24 ,  BAL46 ,  BAL68  and  BAL810  have  been  averaged 
over  all  gridpoints  for  the  four  computational  dates. 
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The  computed  values  of  BAL24  and  BAL68  were  then  compared  with  one- 


half  of  BAL26  and  BAL610  respectively  to  determine  the  relative  per¬ 
centage  of  heating  (cooling)  in  these  layers,  as  depicted  in  Table  X 
for  the  mid-seasonal  winter  and  summer  dates.  It  is  noted  that  the 
values  of  BAL46  and  BAL810  are  identical  in  magnitude  but  of  opposite 
algebraic  sign  to  BAL24  and  BAL68  respectively. 

The  results  of  these  four-layer  computations  (Table  X)  indicates 
that  in  layer  (2,4)  there  is  approximately  thirty  percent  less  cooling 
than  would  have  been  deduced  in  the  simpler  two-layer  model.  This 
implies  less  radiationally  induced  instability  in  the  upper  troposphere 
than  that  originally  estimated.  In  the  layer  (4,6)  there  is  compensat¬ 
ing  increased  radiative  de-stabilization.  This  radiative  difference 
in  the  layer  (4,6)  appears  to  be  primarily  a  result  of  the  placement  of 
clouds  in  both  models.  For  instance,  the  upper  clouds  have  tops  at 
k=4,  and  IR  net  flux  acts  as  the  chief  source  of  increased  de-stabili¬ 
zation  in  layer (4,6)  (and  of  decreased  de-stabilization  in  layer  (2,4)). 
Likewise  the  lower  clouds  have  tops  at  k=8;  so,  it  follows  that  IR  net 
flux  acts  as  the  chief  source  of  increased  de-stabilization  in  layer 
(8,10),  and  finally,  of  reduced  de-stabilization  in  layer  (6,8). 

The  tuning  of  solar  reflectances  has  not  introduced  a  significant 
radiative  influence  on  the  results  of  the  de-stabilization  differences 
just  outlined.  Rather  it  was  the  placement  of  clouds  relative  to 
layer  boundaries  that  affected  primarily  the  IR  net-flux  calculations. 
The  results  summarized  in  Table  X  are  then  to  be  considered  valid  if 
the  cloud-parameterization  model  approaches  climatological  reality. 
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WINTER  RESULTS  SUMMER  RESULTS 

2-Layer  4-Layer  2-Layer  4-Laye: 

Layers  model  model  model  model 
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vi.  the  zonal  distribution  of  radiational 

BALANCE  TERMS  OF  THE  OCEAN-ATMOSPHERE  SYSTEM 


A.  GENERAL  INTRODUCTION:  ZONAL  CROSS-SECTIONS 

The  zonally  distributed  cross-sections  of  the  radiation  contributions 
over  the  ocean-troposphere  system  are  presented  in  Figs.  18,  19,  20  and 
21  for  winter,  spring,  summer  and  autumn  seasons,  respectively.  The 
cross-sections,  show  the  results  after  averaging  over  the  four  meridians 
on  the  mid-seasonal  dates  considered  in  this  study.  The  results  are  dis¬ 
played  in  the  format  of  Fig.  9. 

In  obtaining  zonally  distributed  seasonal  means  of  radiative  heating 
rates,  denoted  Q  (<f> )  ,  at  each  five-degree  multiple  of  latitude  <j)  in  the 
range  20S,...  65N,  values  of  each  radiative  parameter  listed  in  Fig.  9 
were  arithmetically  averaged  over  the  four  meridians.  At  <f>  =  65N  there 
was  only  one  contribution  to  the  zonal  average,  while  in  the  Southern  Hemi¬ 
sphere  latitudes  20S,  15S  and  10S,  only  two  values  (on  X  =  125W  and  X  = 

35W)  of  each  parameter  contributed  to  the  means.  At  5S,  there  were  three 
meridional  sets  of  radiative  parameters  entering  into  Q(<f>).  Otherwise, 
there  were  four  seasonal  values  of  each  radiative  parameter  entering  into 
the  arithmetic  mean  Q  (40  at  each  latitude  (J>  of  the  cross-sections,  Figs. 

18 , • . . , 21 . 

Consequently,  near  the  northern  and  southern  boundaries  of  Figs.  18, 

. ..,21,  the  listed  seasonal  values  may  not  be  as  representative  as  in 
mid-latitudes.  Nevertheless,  the  general  equatorial- to-polar  trend  in 
the  radiative-change  terms  appears  reliable  in  both  hemispheres. 
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18(a).  Zonally-averaged  radiational  cross-section  for  tropical  latitudes.  Refer  to  Fig.  9 
.  All  values  listed  are  daily  averages  in  ly  min--*-  and  are  computed  from  data  for  16  January 
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Figure  18(b).  Zonal ly-averaged  radiational  cross-section  for  higher  latitudes.  Refer  to  Fig.  9  for 
key.  All  values  listed  are  daily  averages  in  ly  min-1  and  are  from  data  for  16  January  1974. 
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re  19(a) „  Zonally-averaged  radiational  cross-section  for  tropical  latitudes.  Refer  to  Fig.  9 
key.  All  values  listed  are  daily  averages  in  ly  min"-*-  and  are  computed  from  data  for  16  April 
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20(a).  Zonally-averaged  radiational  cross-section  for  tropical  latitudes*  Refer  to  Fig. 

.  All  values  listed  are  daily  averages  in  ly  min“l  and  are  computed  from  data  for  16  July 
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Figure  20(b).  Zonally-averaged  radiational  cross-section  for  higher  latitudes.  Refer  to  Fig.  9  for 
key.  All  values  listed  are  daily  averages  in  ly  min-1  and  are  from  data  for  16  July  1974. 
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21(a).  Zonally-averaged  radiational  cross-section  for  tropical  latitudes.  Refer  to  Fig.  9 
•  All  values  listed  are  daily  averages  in  ly  min-1  and  are  computed  from  data  for  16  October 
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Figure  21(b).  Zonally-averaged  radiational  cross-section  for  higher  latitudes.  Refer  to  Fig.  9  for 
key.  All  values  listed  are  daily  averages  in  ly  min-1  and  are  from  data  for  16  October  1973. 


A  zonal  cross-section  depicting  the  mean  annual  radiative  distribu¬ 
tion  has  been  constructed,  Fig.  22,  by  arithmetic  averaging  over  the 
four  mid-seasonal  sets  of  results  (Figs0  18,.. .,21). 

It  has  also  been  convenient  in  earlier  sections  (cf.  Tables  III,  IV, 
VIII)  to  compute  "weighted  averages"  with  respect  to  latitudes  in  the 
Northern  Hemisphere.  Thus  a  definition  of  a  weighted-average  parameter 
which  takes  into  account  the  number  of  observations  k.  available  at  each 

l 

latitude  has  been  given  (Meyers,  1975)  as 


14  4 


Here  Q..  is  the  Q-value  on  meridian  j  at  latitude  6.  and  k.  =  1,...,4  is 
j  l  ii 

the  number  of  meridional  observations  available  for  the  arithmetic 
average  Q((JO.  Note  that  i  =  1,...,14  corresponds  to  the  14  latitudes, 
(j)^  =  0,  . . . ,  65N . 

B.  ANNUAL  RADIATIVE  BALANCE  FOR  THE  EARTH-TROPOSPHERE  SYSTEM 

For  purposes  of  this  summary,  the  term  BALT  (at  k=2)  in  Fig.  9  has 
been  redesignated  for  simplicity  as  R^_  (cf> )  ,  and  annual  values  of  R^  have 
been  plotted  in  Fig.  23  as  a  function  of  latitude.  R  represents  the 
net  radiative- transfer  rate  across  the  tropopause,  that  is  at  the  top  of 
the  troposphere-ocean  system  considered  here.  Similarly,  the  net  flux 
at  the  surface  previously  denoted  BALB  in  the  key,  Fig.  9,  is 
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ire  22(a)*  Zonal-annual  radiational  cross-section  for  tropical  latitudes.  Refer  to  Fig.  9 
key.  All  values  listed  are  daily  averages  in  ly  min~^  and  are  computed  by  arithmetic  averaging 
■  the  four  midseasonal  results  (Figs.  18 ,...,21). 
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e  22(b)*  Zonal-annual  radiational  cross-section  for  higher  latitudes*  Refer  to  Fig*  9  for 
All  values  listed  are  daily  averages  in  ly  min-^  and  are  computed  by  arithmetic  averaging 
the  four  midseasonal  results  (Figs.  18, ***,21). 


redesignated  here  following  Malkus  (1962)  as  R(<f>)  and  its  annual  distri¬ 
bution  has  been  graphed  against  latitude  in  Fig,  23.  Finally,  the  rela¬ 
tionship  between  and  R  in  any  column  is  given  by 

R  =  R  +  R  (6-2) 

t  a 

where  R^  is  the  overall  net  cooling  rate  or  flux-divergence  of  the  tropo¬ 
spheric  column,  e.g.. 


R  =  BAL24  +  BAL46  +  BAL68  +  BAL810  (6-3) 

a 

In  Fig.  23,  values  of  R  ($)  have  also  been  plotted  against  latitude;  how- 

a 

ever,  a  simpler  method  of  obtaining  R  from  (6-2) ,  namely, 

a 

R  =  R  -  R  (6-4) 

Ci  Lr 

has  been  employed.  The  resulting  zonal  annual  distributions  of  each  of 

the  three  parameters,  R  ,  R  and  R  ,  are  shown  as  functions  of  latitude 

t  a 

in  Fig.  23,  where  they  are  superimposed  against  similar  functions  drawn 
from  Sellers  (1965) « 

The  comparison  which  follows  focuses  on  the  annual  distribution  in 
the  Northern  Hemisphere  only.  It  should  also  be  noted  that  Sellers1 
radiative  parameters  were  taken  from  climatology  at  the  top  of  the  mean 
zonal  atmosphere ,  whereas  this  radiative  model  gives  corresponding  results 
for  level  k  =  2,  based  upon  tropospheric  soundings  over  an  oceanic  surface 
only. 

Throughout  the  latitude  range  (0-65N) /  the  radiative  net  flux  R  at 
the  ocean  surface  is  substantially  greater  than  the  climatological  amount 
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Figure  23.  Radiative  net  fluxes  at  the  tropopause  (R  )  at  the  ocean 

surface  (R) ;  and  net-flux  divergence  for  the  tropospheric  column  (R_) 

from  the  zonal  annual  results  (Fig.  22(a,b)).  Solid  lines  indicate' 

values  of  R.  R  and  R  after  Sellers  (1965). 
t  a 
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of  Sellers.  This  is  attributable  to  the  reduced  model-albedo  values  of 

Sec.  IV.  F.  This  in  turn  allowed  more  insolational-absorption  in  the 

ocean  thereby  increasing  R  over  all  latitudes  as  compared  to  Sellers1 

values  (Fig.  23).  The  term  R-R  (which  equals  -R  )  may  be  regarded  as 

t  a 

an  atmospheric  "greenhouse  effect"  in  contributing  to  warming  of  the  sur¬ 
face  both  for  the  model  and  for  Sellers1  climatology.  For  each  latitude. 
Fig.  23  indicates  that  the  "model  greenhouse"  is  in  close  agreement  with 
that  of  Sellers. 

The  Redistributions  of  both  Sellers  and  of  the  present  model  are  in 
good  agreement  across  the  entire  range  of  Fig.  23,  both  representing  cool¬ 
ing  rates  in  the  troposphere  for  all  <J>*  A  minimum  value  of  R^  at  latitude 
20N  appears  on  both  curves  and  is  presumably  due  to  a  local  maximum  of 
IR  net-flux  divergence  F^*  -  (resulting  from  the  minimum  opaque 

cloud  cover)  in  subtropical  latitudes.  From  SON  to  65N,  an  increase  in 
the  model  values  of  R ^  appears  in  Fig.  23  and  seems  to  indicate  a  reduced 
IR  net-flux  loss  associated  with  high  CL-values  (Table  III)  at  these 
northerly  latitudes.  Sellers*  results  indicate  the  opposite  trend,  i.e., 

increasingly  negative  R  -values  from  50-60N,  but  his  climatology  may  not 

a 

include  the  kind  of  detail  that  would  indicate  reasonable  oceanic  cloudi¬ 
ness  over  these  latitudes. 

The  radiative  model  gives  R^-values  which  are  generally  larger,  both 
over  tropical  latitudes  as  well  as  over  the  northern  latitudes,  than  the 
R^-values  of  Sellers.  As  in  the  R((f))  comparison  this  is  a  result  primar¬ 
ily  of  the  reduced  model-albedo.  In  the  mid-latitude  zone  25-40N  of 
Fig.  23,  close  agreement  occurs  between  the  two  R^-curves  presumably  due 
to  such  similar  climatological  effects  as  for  example,  the  mean  polar 
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front.  The  positive- to-negative  R  crossover-point  occurs  very  near  to 
40N  for  both  Redistributions  of  Fig.  23.  The  model  weighted-average  is 
R^_  (cj))  =  0.048  ly  min  ^  if  averaging  is  considered  only  to  65N.  If  nega¬ 
tive  R  annual  values  exist  in  latitudes  (p  65N,  an  ocean- troposphere 
radiative  balance  should  very  nearly  be  established  over  the  Northern 
Hemisphere.  However  data  for  the  model  were  not  analyzed  over  ice-covered 
regions  poleward  of  65N,  where  presumably  higher  surface  albedos  and 
cooler  drier  soundings  would  enable  R  to  assume  larger  negative  values. 

The  radiative-model  values  of  R  ((}))  ,  R  (4> )  and  R  (c|) )  are  in  agreement 

t  a 

with  present  albedo  and  IR  net-flux  observations  from  recent  satellite 
climatology.  It  should  be  noted  that  Sellers1  radiative  results  are 
based  on  climatology  of  London  (1957) ,  and  include  also  data  from  the 
land  areas  of  the  Northern  Hemisphere;  on  the  other  hand  the  radiative 
model  values  have  been  based  only  upon  soundings  over  the  ocean. 

C.  CROSS-SEASONAL  EFFECTS  IN  THE  NORTHERN  HEMISPHERE 
1 .  On  the  Net  Flux  Across  the  Tropopause,  R^ 

The  seasonal  distribution  of  the  tropospheric  net  flux  is  pre¬ 
sented  as  a  function  of  latitude  in  Table  XI (a)  where  the  listed  R^-values 
correspond  to  the  tropopause,  level  k  =  2.  One  may  discern  two  geographic 
regimes  of  annual  waves  in  Table  XI (a);  a  polar  zone  (25-65N)  and  a 
tropical-subtropical  zone  (0-25N) . 

The  cross-seasonal  variation  of  R^  in  the  northern  regime  has 
the  following  property:  R^-values  are  a  minimum  in  winter  increasing  to 
a  maximum  in  summer  followed  by  a  decrease  in  autumn.  In  the  tropical 
zone  (0-10N) ,  R^-values  have  positive  peaks  in  spring  and  autumn.  The 
summer  minimum  of  R  in  0-10N  is  presumably  attributable  to  increased 
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cloudiness  associated  with  the  ITCZ  in  this  latitude  zone.  The  subtropi¬ 
cal  region  15-25N  does  not  have  a  minimum  in  midsummer,  and  may  be  simply 
classified  as  intermediate  in  behavior  between  that  exhibited  by  R  in 
the  polar  and  the  tropical  regimes.  These  geographic  cross-sectional 
effects  (Table  XI  (a) )  tend  to  confirm  the  general  validity  of  the  radia¬ 
tive  model. 

2.  On  the  Sea-surface  Model  Balance,  R 

The  seasonal  distribution  of  the  net  radiative  flux  R  at  the 
earth's  surface  is  listed  as  a  function  of  latitude  in  Table  XI (b) .  Again 
a  polar  zone  (25-65N)  and  a  tropical-subtropical  zone  (0-25N)  are  defined 
for  discussing  cross-seasonal  effects  on  R. 

The  model  results  for  R  are  analogous  to  those  just  specified  for 
R  ,  namely: 

(i)  in  the  polar  latitudes  there  is  a  single  sine  wave  with 
annual  periodicity  and  maximum  in  midsummer; 

(ii)  in  tropical  latitudes  0-10N,  the  double  peaked  character 
with  maximum  in  spring  and  autumn  appears; 

(iii)  there  exists  an  intermediate  subtropical  zone  of  R-values 
which  no  longer  exhibits  a  minimum  in  summer. 


D.  TOP  OF  THE  ATMOSPHERE  COMPARISON  OF  MODEL  WITH  SATELLITE  DATA 
1 .  Net  Radiative  Transfer  Rate  at  Top  of  the  Atmosphere 

The  three  parameters  that  describe  the  radiative  flux  at  level 
k  =  0  are  associated  in  the  following  way: 


\  ~  FT 


(6-5) 


Q 


N' 


the  daily  averaged  incoming  solar  flux. 


can  be  related  to  QAVE  and 
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QREF  (Fig.  9)  by  the  relationship 


Qn  =  Q^E  (1-ALB)  =  QAVE  -  QREF 


(6-6) 


where  QAVE  and  OREF  have  already  been  evaluated  according  to  Eqs .  (5-3) 

and  (5-6)  respectively  and  the  factor  2.00/1.92  has  the  effect  of  deriv¬ 
ing  the  mean  daily  solar  insolation  at  level  k  =  0  before  ozone-oxygen 
absorption.  F  ,  the  IR  net  flux  to  space,  was  previously  discussed  in 
Sec.  III.D.  [Eq.  (3-11)],  and  R^,  the  difference  between  Q^  and  F^, 
represents  the  daily-averaged  radiative  net  flux  at  the  top  of  the 
atmosphere.  In  order  to  more  easily  distinguish  between  model-values 
of  R  and  the  same  parameter  from  satellite  climatology  of  Raschke  et  al. 
(1973),  R^  is  further  identified  as  R^MOD  or  R^RAS,  respectively.  QAVE 
is  assumed  identical  both  for  the  model  and  for  the  treatment  of  the 
Raschke  data  in  this  study. 

2.  Zonally-Averaged  Computations  of  R^ 

With  the  use  of  Eqs.  (6-6)  and  (6-5),  values  of  Q^,  F^  and  R^ 
for  each  of  the  four  mid-seasonal  dates  and  by  five-degree  latitude 
intervals  have  been  computed  both  for  the  model-parameters  and  the  satel¬ 
lite-observed  (RAS)  parameters.  The  results  of  the  computations  both 
seasonally  and  for  the  annual  mean  case  are  presented  in  Table  XII. 

A  comparison  of  R^MOD  with  R^RAS  from  the  zonal-annual  results 
at  the  bottom  of  Table  XII  shows  that  the  crossover  from  positive-to- 
negative  values  of  R^  occurs  near  40N  for  both  model-  and  Raschke-results . 

A  statistical  accuracy  check  of  R^MOD  as  against  R^RAS  data  is 
conducted  using  the  standard  deviations  of  the  difference  parameter 
(R^MOD  -  R^RAS)  over  all  Northern  Hemisphere  latitudes.  A  summary  of 
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the  mean  and  standard  deviations  of  this  parameter  over  the  latitude 
range  of  these  five  cross-sections  (Table  XII)  follows: 


16  Jan. 

16  Apr. 

16  Jul. 

16  Oct. 

Annual 

R^MOD-R^RAS 

.004 

.024 

.019 

.014 

.015  ly  min 

Std.  Dev. 

.028 

.037 

.053 

.039 

.039  ly  min 

The  zonally-averaged  mean  (R^MOD-R^RAS)  is  positive  for  each 
season  as  well  as  for  the  annual  distribution  (whose  mean  is  0.015  ly 
min  .  This  small  positive  difference  was  anticipated  from  considera¬ 
tion  of  Sec.  IV. F.  (Table  IX)  where  the  middle- to-high  latitude  albedos 
were  tuned  slightly  too  small. 

The  overall  standard  deviation  in  (R^MOD-R^RAS)  of  0.039  ly  min  ^ 
is  of  the  same  order  of  magnitude  as  that  of  Q^MOD-Q^RAS .  This  compara¬ 
tively  small  difference  in  the  case  of  model  versus  satellite  climatology 
was  a  considerable  improvement  over  (A,B,C,D)  where  untuned  albedos  were 
used. 


E.  STRATOSPHERIC  MODEL  RADIATIVE  BALANCE 

It  is  clear  that  the  flux-convergence  in  the  stratosphere  is  given 
by 


where  R  is  given  by  Eq.  (6-5).  Note  that  R  is  identical  to  BALT  of 
N  t 

Fig.  9  which  reduces  to 

R  =  OAVE  -  QREF  -  F 2*  .  (6-7) 


110 


Forming  the  difference,  R  -  R  .  leads  to 

N  t 

R  -  R  =  .04  CAVE  -  (F  -  F  *)  .  (6-8) 

N  t  ~  T  2 

The  stratospheric  flux-convergence,  R  -R^,  is  displayed  in  zonal- 
annual  format  in  Table  XIII.  Based  on  the  Northern  Hemisphere  weighted 
results  (bottom  line) ,  it  appears  that  a  positive  net-flux  convergence 
exists  between  0-65N.  However  the  trend  in  R^-R^  at  higher  latitudes 
is  towards  increasingly  negative  radiative  values  in  the  stratosphere 
poleward  of  65N  such  that  an  annual  mean  radiative  balance  may  be  in¬ 
ferred.  For  instance,  assuming  that  in  the  latitude  range  65-90N,  an 
average  zonal-annual  value  of  R^-R^  =  -0.0502  ly  min  ^  exists  at  77. 5N, 
then,  the  cosine-weighted  mean  is  -0.0126  ly  min  ^  in  the  region  pole- 
ward  of  65N.  Thus  the  summed  Northern  Hemisphere  weighted  mean  would  be 
zero  and  a  radiative  balance  would  exist  in  the  stratosphere. 

F.  ZONAL- ANNUAL  NORTHERN  HEMISPHERE  HEAT  BUDGET 
1 .  Tropospheric  Radiation  Budget 

The  radiative  heating  rate  of  a  tropospheric  column  may  be  ex¬ 
pressed  as  a  function  of  R  (<$>)  by  the  right  side  of  Eq.  (6-9) 

a 

Q  +  S  =  R  +  (E  +  Hr)  (6-9) 

va  a  a  i 

is  the  storage  heating  rate  of  the  column  and  Q  is  the  required 
flux-convergence  of  sensible  and  latent  heat  energies  compatible  with 
the  heat  balance  at  latitude  cf>.  R^  is  the  tropospheric  radiative  net 
cooling  rate  and  is  depicted  in  Fig.  23.  (E  +  Hp)  includes  the  latent 
and  sensible  heat  parameters  respectively,  but  these  parameters  are 
not  part  of  this  study  and  will  be  considered  to  be  contained  in  Q 
of  Eq.  (6-9)  for  simplicity. 


Ill 


Lat. 

R 

N 

Rt 

\  " 1 

2  OS 

.1179 

.1030 

.0149 

15 

.1260 

.1085 

.0175 

10 

.1293 

.1074 

.0219 

5 

.1456 

.1248 

.0208 

0 

.1443 

.1215 

.0228 

5 

.1337 

.1097 

.0240 

10 

.1230 

.0969 

.0261 

15 

.1084 

.0826 

.0258 

20 

.0907 

.0644 

.0236 

25 

.0603 

.0367 

.0236 

30 

.0441 

.0293 

.0148 

35 

.0221 

.0113 

.0108 

40 

.0081 

.0052 

o  0029 

45 

-.0137 

-.0107 

-.0030 

50 

-.0273 

-.0190 

-.0083 

55 

-.0426 

-.0168 

-.0258 

60 

-.0566 

-.0290 

-.0256 

65N 

-.0563 

-.0106 

-.0457 

Wt.  Avg. 

.0602 

.0476 

.0126 

TABLE  XIII.  Computed  radiative  net  fluxes  at  the  top  of  the 
atmosphere  (R^)  at  the  tropopause  (R  ) ;  and  the 
stratospheric  absorption  (R  -R^) •  All  parameters 
are  listed  in  the  zonal-annual  format  in  ly  min"1. 
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The  Northern  Hemisphere  cosine-weighted  turns  out  to  be 

R  =  -  0.1551  ly  min  ^ 
a 

and  if  it  is  temporarily  attributable  to  mean  storage-cooling  S  of  the 

a 

troposphere  alone,  the  mean  storage  rate  corresponds  to  a  temperature- 
change  rate  given  by 

6t  1440  <V 

<6t}  =  4.1  -  °C/ day  (6-10) 

mb 

with  AP^  =  800  nib  in  the  troposphere.  The  resultant  cooling  rate  (6-10) 

over  the  tropospheric  depth,  800  mb,  with  zero  lateral  flux  divergence, 

2 

is  approximately  -1.14°C  per  day  averaged  over  the  mean  cm  tropospheric 
column. 

Thus  the  general  circulation  of  the  atmosphere  would  have  to 

function  to  bring  about  the  atmospheric  heat  balance.  This  process 

would  require  the  proper  flux-convergence  Q  to  offset  the  annual  radia- 

2 

tive  loss  of  the  average  cm  column  of  the  troposphere. 

2 .  Zonal-annual  Heat  Budget  of  the  Ocean 

An  equation  similar  to  (6-9)  may  be  written  for  the  zonal-annual 
heat  budget  for  the  ocean  as  a  function  of  R(cj))  . 

Q  +  S  =  R  -  (E  +  Hr)  (6-11) 

vo  o  r 

Here,  S  is  the  storage  rate  of  the  ocean  water  mass  and  Q  is  the 

O  vO 

required  mean  oceanic  heat  flux  divergence  for  a  balance  at  latitude  <J). 

The  primary  ocean  mass  heating  parameter  is  the  net  radiative  heating 

flux  at  the  ocean  surface,  R  (Fig.  23),  and  (E  +  Hp)  is  for  convenience 

considered  part  of  O 

*vo 
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The  Northern  Hemisphere  annual  weighted  mean  value  of  R  is 
R  =  0.2027  ly  min  1 

which  corresponds  (Table  XI)  to  a  mean  radiative  heating  rate  in  the 

water-mass  column.  It  should  be  noted  that  this  mean  heating  rate  is 

noticeably  larger  than  R  and  of  opposite  sign  to  the  corresponding 

a 

tropospheric  radiative  cooling  effect  (R  =  -0.1551  ly  min  0  How- 

a 

ever,  if  the  three-dimensional  oceanic  flux  divergence  of  both 
sensible  and  latent  heat  is  considered,  the  oceanic  heat  balance  should 
result. 
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VII.  CONCLUSIONS 


The  major  change  in  this  radiative  model  from  the  model  evolving 
from  (A,B,C,D)  was  the  introduction  of  tuned  cloud  reflectances.  This 
innovation  brought  the  resultant  model-albedos  into  agreement  with 
satellite  climatology  for  comparable  data  dates. 

Gridpoint  comparisons  at  the  top  of  the  atmosphere  of  model  net  flux 
and  of  the  net  incoming  model  insolation  were  made  with  similar 
parameters  from  Raschkes 1  satellite  climatology  (1973).  Closer  agree¬ 
ment  was  obtained  between  Q  MOD  and  Q  RAS  than  in  any  of  the  preceding 

N  N 

studies  (A,B,C,D) .  Similarly,  close  agreement  in  a  least-squares 
sense  was  achieved  here  between  F  MOD  and  F^RAS  primarily  because  of 
an  improved  formulation  for  the  F  - computation  than  was  used  in 
earlier  studies.  Finally  the  ocean-atmosphere  system  net  flux 

\  “  0N  -  Ft 

has  only  a  small  least-square  error  when  model  and  satellite  results 
are  compared  with  properly  stratified  gridpoint  and  time  cross-sections. 

In  this  study,  mid-latitude  cloud  reflectances  were  also  tuned  to 
obtain  closer  agreement  between  the  global  albedo  and  corresponding 
values  from  satellite  climatology.  The  result  was  to  systematically  in¬ 
crease  relative  to  Raschke's  values.  This  result  was  undoubtedly 
due  to  a  mismatch  between  the  model  clouds  and  those  involved  in  the 
climatological  average.  The  conclusion  reached  is  that  if  tuning  of 
cloud  reflectances  of  this  nature  is  to  be  attempted  in  mid-latitudes 
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using  FNWC  soundings,  more  accurate  cloud-parameterization  formulas 
should  be  established  against  observed  cloud-cover  amounts  for  identical 
experiment-periods  similar  to  GATE .  If  then  further  cloud-reflectance 
tuning  is  necessary,  it  could  proceed  from  a  more  certain  knowledge  of 
the  "ground  truth"  provided  by  these  experiments. 
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APPENDIX  A.  COMPUTER  PROGRAM 


//FORT .SYS  IN 
COMMON 


502 

500 


C 

r 

c 

c 

c 


CD  * 

6T(  11  )  ,DA74( 11 , 5)  ,CL (2  ) 

CC MMO N/ ABS /  FA2 » AL  PH AG  *  C Z , 5  ECZ , F S , FA  D J 
COMMON/ 7  k  I G/  3r i C7,SD,CD, AL  AT 
COMMON/ VP/  Q( 11 )  , QS (11  ) 

COMMON/ ARM/  RM 
COMMON/ R ASK C/RAS  VAL 

CO MMO N/ CASK/  AL  BRAS (93) ,F2RAS( 53 ) , ISO 
COMMON/ FACT/  K,ISEA 
FQRMAT(4X,2(F8.4,2X> ) 

FORMAT  (2X,F6.1,lX,F7.2,lX,5(F9.5,iX),<>X,2I2) 
DO  63  ISEA=1,4 

■Q.l)  DEC=-21. 07917/57. 29578 
'  “  DEC  =  8.43333/  57.29578 

DE  C  =  21.5000/57.29578 

DE  C  =  -8.225/57.29576 

RM  =  . 9837 1 5**2 
RM  =  1 • 30  33  3**2 
R M  =  1  . C  1 644**2 
IF (ISEA. EG. 4 >  RM  =  . 997165**2 


I  F ( I SE  A . 
IF  (  IS  EA. 
IF  (ISErt. 
IF  (  ISEA. 
IF (IS EA  . 
IF  (  ISEA. 
I  F (  IScA 


EQ .2  ) 
Eg. 3 ) 
EQ.4) 
Ey.l ) 
6Q.2  5 
EG.3  ) 


THE  FOLLOWING  READ  : 
CLIMATOLOGY  DATA  Cp 
ROUTINE  REFT. 


IT AT EM ENT 
RASCHKE 


= 


APPLIES  TO  THE  ... 
ET  AL.  AS  USED  IN 


TELL  I  7  F 
7F!c  S  0  B- 


1 


1C 


ISEA. SO.?. AND. ISC. EQ. 25)  L 
L  ANGLE  (C-Z,  SEC  Z,  i  3  0,  LABI  ,  D 
ISEA.EQ.l .AND. ALAT.LE.25. ) 
ISEA.EQ.l . AND. ALAT.G7.25. ) 
ISEA  .  EQ . 2 . AND .ALA l  . L  —  . 25. ) 


_\D.ALA_ 
.  Ar>  D  .  A  L  A  i 
AND. ALA  I 
EQ.^5  .  AND  .Ala  i 
EQ.4. AND.ALAT 
EQ.4. AND. ALAT 


.  EQ  •  2 
,  EQ .2 

i  c  Q  . 


Li 
G  t 

,LE, 
i  GT 
LE, 
.  GT 


25 
25 .  ) 
25.  ) 
25.) 
25.  ) 
25.  ) 


_ =  24 

EC) 

RAS VAL 
RASVAL 
RASVAL 
RASVAL 
RASVAL 
RASVAL 
RASVAL 
RASVAL 


,  H 
,8 
i  3  5 
,65 
,  6 
,  7 
,  4 
,6 


20 

66 

68 

95 


CC  9  130  =  1 » 92 

RE  AC (  5,  502)  AL6R A S ( I  SO ) , F2R AS ( I  SO ) 

CONTINUE 
CC  66  150=1,93 
DO  10  1=1,11 

F SAD( 17  ,500, =ND=99  )  ( DATA ( I ,K),K=1,2),0(I),GS(I) 
(DATA (I ,KK) ,KK=3, 5) ,LAB1,LAB2 
CCNT  INUE 

IF  (ISEA.EQ.3.ANQ.ISC.EQ.25)  LA- 

CALL  ANC'  '  -  J  ‘ 

IF  (  ISEA 
IF  (ISEA 
IF (ISEA 
IF  (ISEA , 

IF  (  ISEA 
I F ( IS  E  A , 

IF (ISEA, 

IF (ISEA 
C  <LL  BB E 

ALPFAG=(.54*(.7-CZ) )  +  •  06 
I  F ( ALPHAS. LT..C6 )  ALD|-AG=.06 
FS=.65i*CZ*1.92/RM 
F  ACJ=  (2.00-CZ  )/RH 
F  A  2= ( .345*1. 92*CZ ) /RM 
CO  20  1=1,4 
CL ( 1) =0.0 
CL ( 2 ) =0  .0 

IFU.EQ.2.0R.I.BQ.4)  CL(1)  =  1.000 
IF{I.GE.3  )  C L ( 2 )  =  1 . 00 0 
WT=(1.0-DATA(6,3) )*( 1 .O-DATA ( 2 , 2 ) ) 

I  F (  I . E Q . 2 )  WT  =  DA 7  A ( 6 , 3 ) * ( 1 . O-CA  T  A ( 2 , 

I  F ( I • EQ  .3 )  WT=DA7 A  (  2, 3 ) ^  '  ’  "  r  ’ 7 '  " 

IF (I. EQ.4)  WT  =  DAT  A ( 2 , 3 ) 

CALL  A6SCJRB(  CL  ( 1  )  ,CL(2  ) 

CONTINUE 
CONTINUE 
CONTINUE 
STOP 
END 

SUBROUTINE  A3S0RB ( Cl , C2 , WT ) 

CCMtMON/ABS/  FA2,  AL  PHAG,  C  Z ,  S  ECZ  ,  FS  ,  FAD  J 
C.  C  M MON  /  R  A  S  KC  /  R  A  S  V  A  L 
1  F ( Cl  .GT. .5.0R.C2.GT..5)  GO  TC  20 
ALPHAft= ( .085+( .25074*ALOG10 (SECZ ) ) ) 

IF (ALPHAR.GT. 1.000)  ALPHAR= 1.000 


/  '"•  l  1  •  W  —  UH  I  Hid,  j)  ) 

)*(1.G-D4TA(6»3) ) 
*  *DAT A ( 6 , 3  ) 

,  WT  ) 


117 


A  = (  ( 1 • - A  L  PHAR ) * ( 1 .-ALPHA G) ) / d.- <  ALPHAR*  ALPFAG  ) ) 
SOL=-A*F  S 
GC  TO  21 

2C  R1=.54*RASVAL 
R  2  =  • 66*RA  S VA  L 
IP(C1.LT..5)  R 1  =  0  •  0 
I F ( C2 . LT. • 5)  R2-G.G 
A= ( 1.-R1 )*(1 .-R2 >*( 1.- ALP HAG) 

C=1.0--(  (Rl-"‘c2)  +  (R2*ALPHAG)+(Ri*ALPHAG)  ) 

C  =  2  •  0  *  R 1  *  R  2*  A  L  P  h  A  G 
SCL=( A  *  F  S ) / ( D+C ) 

21  !F(Cl.L7..1.AND.C2.LT..l)  CALL  AEOO< WT, SCL  ) 

I F ( Cl.GT. . 9.AND.C2. LT. . 1 )  CALL  A61Q(WT,SCL) 

IF  (Cl.LT..  1.AM0.C2.GT..9)  CALL  A  BO  1  (  WT  ,  SCI.  ) 
JF(Ci.GT..9.AND.C2.G7..9)  CALL  AB11(W7,SCL) 
RETURN 
END 

SUBROUTINE  TD<  ANS, IND,N1,N2 ,SEC, C) 

CCNMON  BT(  11  )  ,DATAdl,  5),CL<2) 

U  =  DATA(M,4)-  DAT  A ( N2  »  4 ) 

A =.203 

ANS=C*(  (U*SEC  >**A  ) 

I F ( IND.G7.10. )  AN  S  =  1 . O-ANS 

RETURN 

Ei\C 

5 L GROUT  INF  BBB 

COMMON  BT( II ) , DAT  A  d l »  5 ) ,CL ( 2  ) 

£  I GMA  =  1 .170403 
E= 1 • E-7 
CO  10  1=1,11 
A=CATA( I ,2 )+273.16 
BT ( I }=( A**4) *SIGMA*E 
10  CONTINUE 

C  AT  A  (  2  ?  4  )  =  ( DAT  A.  (  1 , 4  )  +OATA  (  3 ,4  )  )  /2 .00 
D  A  T  A ( 2 , 5 )  =  ( DAT  A ( 1 , 5 ) +DATA ( 3 , 5 ) )/2.CC 
RETURN 
END 

SUBROUTINE  ANGLE (CZ,SECZ,IS,L2,D) 

CCPKGN/TRI 0/  £7, CT, SO, CD, AL AT 

R=57. 29578 

h=35.0/R 

IF(I£.LE.25)  H=55.0/R 
I F ( IS  .GE.26.AND. IS.LE. 50)  H  =  1C.0/R 
A=973.752 
S=  2.0 

I F( IS. GE. 26. AND. IS.LE. 50)  S=I. 0 
AK=L2+1 

IFdS.GE.68)  AK=63-L2 
B={S*((32.-AK>**2)  ) 

$7= ( A-6 ) / ( A+S ) 

CT=(1.- (ST**2> )**.5 
SC=SIN(C  ) 

CC=COS( C) 

CF=COS ( H ) 

CZ=(ST*SD)+(CT*CD*CH) 

SECZ=1.00/CZ 

TAM=ST/CT 

ALAT=( ATAN(TANT) )*57. 29578 

RETURN 

ENC 

SUBROUTINE  ABOO(WT,SOL) 

CGMMGN/A6S/FA2, ALPHA G,CZ, SEC Z, F$, FAD J 
CALL  TD(U24,0,9, 7,SECZ,.271  ) 

CALL  TC(U26,0,9,5,SECZ,.2.71) 

CALL  TC(U23,0,9,3,$5CZ,.27i) 

CALL  TC(U21C,0,9,1,SECZ,.271) 

CALL  -TC(U10,12,9,l,StCZ,.271) 

A24=FA2*U24 
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A46=FA2*<  U26-U24  ) 

A68=FA24(U28-U26  ) 

A610=FA2*(U210~U28> 

T  R ANA  =  F A  2*U10 
A  I  1 0  =  T  R  ANA*<  1.- ALP  HAG  ) 

R EFA=FA2-A24~A46-A6S-A81C~AI 10 

CALL  REFT ( 1, A24, 446, A68, A81C, A  I  1C, TRANA , REFA ,S0L , WT  ) 

RETURN 

END 

SUBROUTINE  A  810 ( WT , SOL ) 
CCMM3N/ABS/FA2,ALPHAG,CZ,SECZ, F$  ,  FAD  J 
CCMMCN/ RASKC/R ASVAL 
R 1  =  • 464  R  AS  VAL 
A 1=. 20 

CALL  TC(TD24,12,9,7,SECZ,.271) 

CALL  TD(TDo8,12, 5, 2,1  ,6667,  .271  ) 

CALL  TO (TD61 0,1 2, 5,1, 1.66667, .271) 

CALL  TD (T 08 10,12 ,3,1, 1.6667, .  271  ) 

F4D=FA2*TD24 
F40=R 1- F4D 
F2L  =  F4li*TD24 
A24=FA2-F40+F4U-F2U 
A46=F4D*A1 
F6C=F4D*<  1.-R1-A1  ) 

F8C=F6C*T068 

F10D=F6D*TD610 

F10U=Fi CD4ALPHAG 

F6U=F1CU*TDS1G 

F6U=F10U*T0610 

F6CD=FfcU*Rl 

FfiCC=F6C0*TD68 

F 10DC=  F6DQ4TD610 

Ad£  =  F6D-FSD+P8’J~F6U+F6DD-F80D 

A £10=  F8C-F10D+F10U-F8U+F8DS-F10D0 

AI10=(F10D4{ 1.- ALPHAS) )+F10DD 

TRANA=F 1CD+F10DD 

R6FA=FA2- A24-A46- A68-A8 10-A ! 10 

CALL  REFT ( 2, A24, A46, A68 , A8 10 , A  I  1C , TR ANA , RE F A , SOL , WT ) 

RET  URN 

ENC 

SUBROUTINE  A301(WT,S0L) 

CCMMCN/ABS/FA2, ALPHA G,CZ, SECZ, FS , FAD J 
COMMON/ R AS KC/R AS VAL 
R  2= . 504 R A  S  VA  L 
A2=.30 

CALL  TD(TD24.12,9, 7,SFCZt.271 ) 

CALL  TC(TD26,12,9,5,SECZ,.271) 

CALL  TC(TD28,i2,9,3,SECZ, .271) 

CALL  TD( 7  0910,12 ,2,1, 1.6667, .271  ) 

CALL  TC(TD68,12, 5, 3, SECZ, ,271) 

CALL  Tu( TD48*i2,7, 3, SECZ » .271 ) 

FAD=FA24T024 

F6C  =  F A  24TD26 

FcC=FA2*7D23 

F  EU  =  F  8C  *R2 

FfcL=F8U4T068 

F4L=F8U*TD48 

F2l=F8U*TD28 

A24  =  FA2— t-40+F4U— F2U 

A4fc=F4C- F6D+F6U-F4U 

A68  =  F6C-F8D«-F8U-F6U 

F  9C  =  F  804 ( 1.-R2-A2) 

F1GD=F9C*TD910 
F10U=F1CD*ALPHAG 
F9L=F1GU*TD910 
F  9  C  D=  F  9  U  *  R  2 
F10CC=F9DD*TD910 

A£10=( F8D4A2 ) +F9D-F 10 D+F 10 U-F9U+ F 9DD-F1 0 DD 

TRANA=F 10D+F10DD 

AI10=(F10D*( l.-ALPHAG) )+F10DD 
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P  E  FA=  FA  2-  A 2 4- A 4 o- A  6 6~  A  9  10- A  I  10 

CALL  Rc  FT ( 3 , A24 , A46 , 46  6 , A8 10, A  I  10 , TRAN  A , RE  FA , SOL , WT ) 

RETURN 

END 

SUBROUTINE  ABii(W7,SCL) 

C C  M.v.ON  /  A  8  S /  F A2  .  A  L  P H A G ,  C  Z ,  S EC Z  t  F S  ,  F AD  J 

CCMMCN/RaSKC/RASVAL 

R1=.46*RASVAL 

R2«.50*RASVAL 

A  1  =  • 20 

A2=.30 

CALL  TC(TD24,12f 9»7,SSCZ».27i) 

CALL  T C ( T 068 ,12, 5,3,1.6667,  .271 ) 

CALL  TCdDSlO,  1 2 , 2 , 1 ,  1 .  o66  7  ,  *  27  1  ) 

F4C  =  F  A  2^~  D  24 

F4U  =  R 1*  F4D 

F2U=F4U*T024 

A24=FA2-F40+F4U-F2U 

A  4  6=  F40*A 1 

F  6C=  F  4D  7-  (  1.-R1-A1  ) 

F8C=F6G*TD68 

FEL=F8C*R2 

F6l=F8U*TD68 

F6CD=F6U*R1 

F8DD=F6DD*TD68 

A68=F6C-F8D+F8U-F6U+F6DC-F3DD 
F9C=F8D*( l.-R2-A2i+F8DD*< 1.-A2  ) 

F10D=F9C*TD910 
F10U=F10D#ALPHAG 
F  9  L  =  F 10  L  *  T  D9 1 0 
F9DD=F9U*R2 
F 1CDC=F9DQ*7D91G 

A8 10= ( ( F  8D  +  F800 ) *A2 1+  F9D-F1CD+F 1CU-F9U+F9DD-F 10DD 
Ai 10= ( F!OQ*< l.-ALPHAG) )+FluDD 
TPANa=F10D+F10DD 
REFA*FA2-A24-A46-A68-A810-AI10 

CALL  kEFT(4,A24, A46,A68,A810, A I  10, TRAN A, RE  FA , SOL , WT ) 

RETURN 

END 

Ste POUT  INF  REFT ( NSUB, A24 , A46 , A68 , A 8 10, A  I  10 , TRAN A , 
1REFA  ,SOL , w7  ) 

COMMON  ET ( 11 ), DATA (11, 5) ,CL(2) 

COMMON/ABS/  FA2, ALPHAG,CZ, SECZ, FS, FAOJ 
CCMMON/TR I G/  ST , CT , SO , CD , AL AT 
CCMMCN/ ARM/  RM 

COMMON/ RA3K/  ALBRAS  (  93  )  ,F2RAS(93)  ,ISO 

CCMMC N/PASKC/  RASVAL 

CCMMCN/ FACT/  K,ISEA 

I F ( NSUB . GT • 1 )  G3  TO  10 

C  A  24=0 • 0 

CA46=0. 0 

C  A63=  0 . 0 

CAE  10=0 . 0 

CAI 10=0.0 

CS  110=0.0 

C  R  E  F  =  0 . 0 

10  CA24=C/,24+A24*tfT 
C A  46  =  C  A  4  6  +  A46*WT 
CA68=CA68+A6S*WT 
CA£10=CA810+A810*WT 
CAI 10=CA1 10+AI10*WT 
CSI 10=CSI 10+SDL*WT 
CR  EF  =  CR  E  F+ ( ( FS-SOL ) *w T ) +R EF A*WT 

CALL  IR (NSUB, A10, A8, A6, A4,A2, FID, F8, F6, F4, F2, WT, FT  ) 

I F ( NSUB • L  T . 4 )  RETURN 
AEG=(CAI10+CSI101/FADJ 
T  AND=  SC/CD 
TANT=ST/CT 
T  AN  =  -TAND*TANT 
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H=ARCOS (TAN) 

S l NH=  S  I  N  (  H  ) 

BARCCS=  (  (  ST  ^-SD#  H  )  +  (C  D*CT  *S INH )  )/3. 14159 
F  200= 1 . 92/  RM 
0  A  VE= F200*3ARC  CS 
A  F  2  =- F  2 

CF  5F=  (CR5F*QAVE) /  (  FAD J*. 9600  ) 

B  ALT  =  CAVE-QREF-F2 
0  24=  <  CA  24*QA  VE  )  /  (  F  AD  J*.  9600  ) 

A  F24= -( F2-F4 ) 

BAL24=C24+AF24 

Q46=  ( CA  46*Q AVE ) / ( FAD J*. 9600 ) 

A  F 4 6= - ( F4-F6 ) 

BAL46=G46+AF46 

06  6  =  ( CA  6  8*  3A  VE ) /{ FAD J*.96C0 ) 

A  F6  S  =  - ( F6-F8 ) 

B  A  L  6  8=  Q68  +  AF68 

0810= (CA81C*QAVS  >/( FAD J*. 9600) 

AF810=-(F8-F10) 

8  AL  810=  08  10  +  AF810 
CA  8G=  { ABG *QAVE) /.9600 
B  A L  6=  CABG-Flu 
AF10=-F10 
CA26=CA24+CA46 

026  =  (  ,5*CA26*GAVE)  /( FAD J*. 9600 ) 

A  F  26=- . 5*  ( F2-F6) 

GAl26=C2o+AF26 

CA61C=CA68+CAB10 

Q6 10= ( . 5*CA610*QAVE)/(FADJ*.9600 ) 

AF610=-.5*(F6-F1Q) 

6AL610=C61C+AF610 

CIF24=B^L24-8AL2o 

DI  F46=BAL46-SAL26 

DI  F66  =  BAL68-BAL610 

CIF810=BAL310-EAL610 

AL  E  =  CRt  F/FAD J 

ALBDI F= A LB-ALB RAS ( ISO) 

F  2  C I  F  =  FT  - F  2R  AS (ISO) 

RNRA$=(CAVE*1»G41667*(1.-ALBFA$ ( I  SO ) ) ) -F2R AS ( IS  0 ) 

R  N  iv  0  C  =  <  CAVE*1. 041667*  (  1.-AL0)  )-FT 

D  I  FRN  =  P  NMCD-RNRAS 

RETURN 

END 

SUBROUTINE  IR(NSUB,FlG,F8,F6,F4,F2,SlG,c8,E6,£4,E2, 
1UT  ,ET ) 

IF(NSUB.GT.l »  GO  TO  10 
OF  F  =  1 .0 
810=0.0 
EE=C.O 
E  6=0 .0 
E4=0.0 
1 2=0 • 0 
ST =0*0 

10  CALL  FF101 Flo, OFF ) 

CALL  FF  8 ( F8  ?  OF  F I 
CALL  F F 6 ( F 6 »  0 F F ) 

CALL  FF 4 { F 4 »  OF F ) 

CALL  FF2( F2»  JFF ) 

CALL  FFTOP ( FT  »  OF  F ) 

E10=E10+WT*F10 

E8=68+WT*F9 

E6=E6+wT*F6 

E4  =  E4tWT  *F4 

E2=E2+WT*F2 

E7=ET+WT*FT 

I  F (NSUb.LT. 4)  RETURN 

E1G=E10/1 440 . 

E6=E8/144G. 

E6=E6/ 1 440 • 

E4  =  E4/ 1440 • 
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E2=E2/1440. 

ET=E 7/1440. 

RETURN 

ENC 

C  SLEROUTINE  FF 10 { F10, OFF ) 

COMMON  BT(  11 )  ,  DATA(  11,5)  ,CL{  21 
I  F  ( CL  ( 1  ).GT..5.0R.Cl<2).GT..5)  C-0  TO  10 
CALL  EWC < EW8,GFF , 3, 1 ) 

CALL  ENC( £W6,GFF, 5,1) 

CALL  EW  C  (  E  W4 , 0  F  F  ,  7 , 1  ) 

CALL  EWC  (  EW2  ,'JFF  ,9,1  ) 

CALL  EWC (EW1, OFF, 10,1) 

CALL  EWC( EW9,GFF, 2,1  ) 

CALL  BOUND ( WAVE , 0  FF , 1 1 , 1 ) 

0 1= ( BT ( 3  )  -  37  (  5  ) )*(EW8+EWo) 

B  2= ( BT (  1  )  —  BT ( 3 ) ) *EW3 
A  1  =  (  (  BT  (  1 )  -87  (  5)  ) - . 5  * ( 31+92) ) 

B1=BT( 10}*WAVF 
B  2  =  ( BT (  9 ) - BT ( 1 0 )  )*(EW2+EW1) 

E  3= ( BT  (  7 ) - BT ( 9 ) ) *  (  £ W4  +  EW 2 ) 

B4=(BT(5)-3T(7) )*<EW6  +EW4) 

A 2=  (BT(5) -(.5* (81+32+33+64) )  ) 

Bl=(i.-(.5*EW9) ) 

A  3  =  ( BT ( 1 ) - BT ( 2 ) )*B1 

1C  F1Q*U l.-CL<  2) )*A1 )  +  ( ( l.-CL<  2) ) * (l.-CL( 1) )**2)  + 

1 (CL  (2)*A3  ) 

RETURN 

ENT 

r 

SUBROUTINE  FF8(F3,QPF) 

COMMON  ET(11),DATA(11,5)  ,  C  L  (  2  ) 

I F  <  CL  < 1 ).GT..5.0R.CL (2) .GT..5)  GC  TO  10 
CALL  EWC ( E  W68 , OF  F , 5 , 3 ) 

CALL  EWC(EW48,CFF,7,3) 

CALL  5 W C ( E W2 8 , OF F , 9 , 3  ) 

CALL  EWC(EW18, OFF, 10,3) 

CALL  EWC ( cW810,OFF, 3, 1  ) 

CALL  BCUND ( WAVE , OFF , 1 1 , 3 ) 

B1  =  EW63^(BT( 3)-3T( 5)  ) 

B  2  =  (  F  W  6  8  +  E  W4  8  )  *  ( 37 (5) -37 ( 7)  ) 

B3=(EW48  +  EW2  8)*( BT( 7) -BTC  9)  ) 

B4=(EW28+5W18)*(3T( 9)-BT<10) ) 

85=  WAV  E*3T ( 10 ) 

A 1  =  6T (3 )-( .5*(61  +  6  2+33  +  B4  +  B 5)  ) 

A2=(BT< 1 ) - BT  C  3 ) ) *( 1.- ( . 5*EW81C)  ) 

A3=(BT( 3)-BT( 5) )*( l.-< .5*EW66)  ) 

10  F8  =  ( ( 1 . -CL ( 1  )  )*A1 )  +  ( (l.-CL(l) )*<1.-CL(2I  )*A2)+<CL(1J 
1*A2)+(CL(  1)*(  l.-CL  (2)  )*A2) 

RETURN 

ENC 

C 

SLBRCUTINE  FF6(F6,0FF) 

COMMON  BT  C 11 J , DATA (11, 5) ,CL(2) 
IF(CL(1).G7..5.0R.CL(2).GT..5)  GO  TO  10 
CALL  6  W  C ( E  W6  8 , OF  F , 5 , 3  ) 

CALL  EaC( EW46,0FF , 7, 5 ) 

CALL  EWC(£W26,0FF,9,5) 

CALL  EWC(EW16f0FF,10,5) 

CALL  EWC ( EW6 10 , 0 F F  ,  5 , 1  ) 

CALL  3CUND(WAVE, OFF, 11,5) 

B 1  =  W  A  V  t  ^  8  T (10) 

G2=(3T( 9 ) - B  V ( 10)  ) * < EW26+EW16) 

B3=(3T( 7) -BT (9 ) )*( SW46+EW26) 

34=EW46*(BT(5)-BT(7) ) 

85=EW68*(BT(3)-BT( 5) ) 

A1=(6T(3)-(.5*(B1+B2+B3+B4+B5) ) ) 

A2= ( BT(  1 )-8T( 3 | ) * ( 1 . - ( . 5* ( 5 W63  +  EW610 ) ) ) 

A  3  =  ( 1 • - ( . 5*EW68)  ) * ( BT ( 3 ) - 3T ( 5 )  ) 

10  F6=(  (  l.-CL  ( J. )  )*A1  )+(  ( 1  .-CL(  i  )  )  *(  l.-CL(  2)  )*A2)+(CL(  1 ) 
1*A3H-(CL(1)*(1.-CL(2)  )*A2) 

RETURN 
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ENC 

C 

SUBROUTINE  FF4  (F4,0FF) 

COMMON  BT( 11) ,DATA( 11 , 5) ,CL (2) 

IF(CL(1  ).GT..5.0R.CL(2)  .GT..5)  GO  TO  10 
CALL  EWC(EW4b,OFF,7,i5  ) 

CALL  EWC( EW48,0FF , 7, 3 ) 

CALL  EwC< SW24,0FF, 9,7 J 
CALL  EWC(EW14, OFF, 10,7) 

CALL  EWC(£W410,0FF,7. 1  ) 

CALL  8GUN0(WAVE,0FF,11,7) 

WAV6=WAVE*57 ( 10) 

61*(BT( 5  )-6T ( 7 )  i -EW46 
E2=(BT  (  3  )  -  8  T  (  5  )  ) * ( EW46+EW48 > 
fi  3= ( BT ( 7 ) - BT ( 9 )  ) *EW24 
B4=(BT (9J-BT (10) ) *( EW14+EW24) 

A1  =  BT  (3  )-(  ,5*(  81 +82+83 +B4+WAVE  )  ) 

A2  =  6T ( 7 )- ( . 5*( 83+ 8 4+ WAVE) ) 

A3=(bT(  l)-BT(B) J  *( 1„-U5*(EW48  +  EW410)  I  ) 

10  F4=( < l.-CL (1 ) )*A1 )  +  (CL(l)*A2 )  +  (  < l.-CL( 1 )  )*{ l.-CU  2) ) 
1*A  2  ) 

RETURN 

ENC 

C 

SUBROUTINE  FF2(F2,0FF) 

COMMON  BT  1  11 ) ,l)A7A( 11,5) , CL<2 ) 
IF(CL<li.GT..5.C'R.CL(2),GT..5)  GO  TO  10 
CALL  EWC ( EW2*,GFF, 9, 7 ) 

CALL  EWC  (  EW26,f)FF,9,  5  ) 

CALL  EWC(cW26,  OFF, 9, 3) 

CALL  EWC(EWi2, OFF, 10,9) 

CALL  EWC( EW210 ,OFF, 9, 1 ) 

CALL  8  C  U  N  D  (  WA  V  t ,  0  F  F ,  1 1 , 9 ) 

C1  =  (BT(  1  >  —  3T  ( 3  )  ) 

C2*(l.-(  (  EW28+EW21O/2.0)  ) 

A3=C1*C2 

61  =  V.AVE-8T  ( 10) 

82=EW12*( BT( 9)-BT( 101 ) 

6  3= ( £W26+SW28)*( BT(3)-3T( 5)  ) 

B4  =  ( E W2 4+ c fa 2 6 i * ( 3T (3)-BT (7 )  ) 

B5=(EW24)£(6T(7)-BT(9) ) 

B=61+ 32+83+84 +8 5 
Al=<67<3)-(.5*8)  ) 

B3--EW24*!  BT (  7)-8T<  9  )  ) 

A  2  = ( 87 (7)-(.5"(Bl  +  S2  +  B3)  )  ) 

10  F2=<A3*<1.-CL<2) )*<1.-CL< 1) ) )+<Al*<l.-CL (1)I)+(A2# 
1CL( 1)  ) 

RETURN 

END 

C 

SUBROUTINE  FFTOP ( FT , OFF ) 

COMMON  BT(il)  , DAT  A ( 1 1 , 5 ) ,CL<2) 

IF(CL(1).G7. .5. OR. CL(2).G7..5)  GO  TO  10 
CALL  EWC(EW01,JFF, 11, 10) 

CALL  EWC(Efc'C2, OFF, 11,5) 

CALL  EWC(EW04, OFF, 11,7) 

CALL  EWC ( E wOb , OFF , 11 , 5  ) 

CALL  EWC(EW03, OFF, 11,3) 

CALL  6WC(EW010, OFF, 11,1) 

B1=<EW01+EW02)*< ST (9) -81 ( 10  )  ) 

B2=( EW02+5W04)*( 3T( 7) -5T(9)  ) 
B3=(EW04+EW06)*(BT(5>-87(7)  ) 

B4  =  <  EW06+SW08)*( 9T<3) -87(5)  ) 

A1  =  (BT (3)-(.5*(Bl  +  82  +  63  +  B4)  )) 

B3=( l.-( .5*( EW08+EWG10 ) ) ) 

A2  =  B3*i BT { 1 1 -BT ( 3  I  I 
A2  =  (  BT  ( 7 ) - ( • 5* ( B 1 +82 )  )  ) 

10  FT  =  <A1*  ( l.-CL(  1 )  )  )  +(  A2*(l.  -CL  (  1  )  )  *< l.-CL  (  2  )  )  )  +(  A3* 
1CL( 1) ) 

RETURN 

END 
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SUBROUTINE  BOUND  (  WAVE  ,CFF,N1, N2 ) 
COMMON  BT( lli ,DATA(11,5) ,CL(2) 

7  =  CAT  A (  10, 2) +273. 16 
U  =  C  AT  A ( N 1 » A ) - 0  A  f  A { N  2 , 4 ) 

IFtU.LT. .00001 )  U= . J0C005 
AL=AL0G10 (U) 

D= . 353* AL-.44 
A  1  =  6. 34  V) **Q 
D=~. 03455* AL-. 705 
A2=U**D 

A3=(6.QG/( .35  3* AL +  3.56)  ) 
W7VE=A1*A2~A3 
IF(OFF.GT.10.)  RETURN 
C  =  CATA(N1, 5)-DAT  A  (N2» 5) 

D  = (  U+.028o) **.26 
81=.07262*(l.-(. 62556*0) ) 

82  =  AL0G10 (C )+( 1 .064) 
h  A  V  E  =  w  A  V  t  +  3 1  *6  2 
RETURN 
END 

SUBROUTINE  E  WC  (  E  W , OFF , M 1 , M 2 ) 
COMMON  31(11), DATA ( 1 1 , 5 ) , CL ( 2 ) 

U  =  DAT  A (  M  1 , 4  )  -  0  AT  A  {  M2 , 4  ) 

D  =  ALC'G1C(U  +  .Q1C) 

A  1= ( • 240*D )  + . 622 
Efc=Al 

IFC3FF.GT.10.)  RETURN 
C  =  CA7  A( Ml , 5) -DAT  A( M2, 5 ) 
81=l.-(.62556*( (U+.0286)**.26) ) 
B2=AL0G1 0(0+1. 0  64 
EU=Al+( .07262*31*62 ) 

RETURN 

ENC 


CONTROL  CARD  FORMAT  FOR  2/3  CLOUD  DATA  AND  "RACIATIVc 
SOUNDING"  DATA  ON  NPS-304 

//GO  .  FT  17F001  DC  UN  IT= 3400-4, VOL  =  $SD =N PS204 , D I $p= (OLD,  K 
//  DCB  =  ( DEN* 2, RECFM  =  FB, LRECL=  77, BLKS I ZE  =  84 H,  LABEL*  (  I,S 
//  DSN=UCGSE7 
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